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Abstract 

The  timing  properties  of  BCF-12  scintillating  plastic  fibers  arranged  in  a  coherent 
bundle  are  investigated.  An  MCNP  transport  model  is  developed  analyzing  the  inter¬ 
action  of  neutrons  with  fiber  bundles  of  varying  size.  A  gradient  in  energy  deposition 
is  seen  starting  at  the  10x10  fiber  bundles  onward.  A  40x40  bundle  has  6.7  ±  0.2% 
intrinsic  efficiency  for  double  scatter  events.  Over  40%  of  the  double  scatter  events 
occur  with  a  time  separtion  of  1  ns  or  greater  in  the  40x40  bundle.  The  lowest  dark 
count  rate  was  achieved  with  the  lens  cap  /  aluminum  can  and  collar  system  with 
a  dark  count  rate  of  15  per  second.  A  bimodal  distribution  of  scintillation  pulse 
decay  constants  is  discovered;  these  two  distributions  are  statistically  different  from 
each  other  allowing  real-time  discrimination  of  dark  events  from  scintillation  events 
using  a  single  PMT.  A  dual-ended  readout  technique  is  used  to  investigate  scintilla¬ 
tion  pulse  shapes.  No  significant  difference  is  measured  between  the  multiple  peaks 
seen  in  scintillation  pulses  of  single  fibers  and  fiber  bundles.  Multiple  scatter  events 
were  not  detected  using  a  7x3  fiber  bundle.  Construction  of  a  40x40  fiber  bundle 
is  recommended  for  future  research  due  to  its  intrinsic  efficiency  and  measurable 
double-scatter  event  probability. 
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INVESTIGATION  OF  BCF-12  PLASTIC  SCINTILLATING  COHERENT  FIBER 


BUNDLE  TIMING  PROPERTIES 

I.  Introduction 


1.1  Background 

Fundamental  problems  with  detecting  large  quantities  of  special  nuclear  material 
(SNM)  at  ranges  of  only  a  few  kilometers  remain  a  national  security  concern.  Nu¬ 
clear  proliferation  and  the  threat  of  nuclear  terrorism  continue  to  grow  at  alarming 
rates.  Consequently,  the  need  for  monitoring  of  illicit  nuclear  trafficking  has  become 
increasingly  important.  Detection  of  SNM  neutron  radiation  at  standoff  ranges  (> 
1  km)  requires  the  capability  to  both  spectroscopically  identify  and  estimate  its  spa¬ 
tial  location  based  on  the  observation  of  very  few  neutron  scattering  events.  Present 
methods  rely  on  cumbersome  and  inefficient  technologies  which  make  them  impos¬ 
sible  to  implement  on  remote  sensing  platforms  such  as  unmanned  aerial  vehicles 
(UAVs)  which  could  patrol  border  crossings  and  shipping  lanes.  A  reasonable  remote 
detector  must  be  reliable,  compact,  efficient,  and  be  able  to  discern  and  extract  signal 
information  from  a  heterogeneous  radiation  environment.  By  using  a  bundle  of  thin 
fiber  scintillators  optically  coupled  to  both  photomultiplier  tubes  (PMTs)  and  charge 
coupled  device  (CCD)  cameras,  an  efficient  neutron  detector  with  inherent  particle 
tracking  spatial  resolutions  can  be  attained.  The  use  of  thin  fibers  allow  neutrons 
to  scatter,  depositing  energy  and  escaping  one  fiber  potentially  depositing  energy  in 
another  fiber,  allowing  for  tracking  of  incident  neutron  energy  and  direction. 
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1.2  Objectives 


The  objective  of  this  research  is  to  investigate  the  timing  properties  of  BCF-12 
scintillating  plastic  fibers  arranged  in  a  coherent  bundle.  Investigating  these  bundles 
is  a  critical  first  step  in  the  development  of  a  novel  time-of-flight  neutron  detector. 
BCF-12  fibers  are  studied  because  of  their  fast  time  response,  peak  wavelength  of 
emission,  and  their  sensitivity  to  neutrons.  These  fibers  are  optically  coupled  to 
PMTs  at  both  ends  for  a  dual-ended  readout  technique  using  digital  pulse  processing 
to  investigate  the  scintillation  pulse  shapes,  spatial  resolution  of  the  fibers,  and  the 
suppression  of  noise.  The  interactions  of  neutrons  with  fiber  bundle  geometries  are 
simulated  to  investigate  the  physics  of  neutron  scattering  in  multiple  fibers.  Results 
from  simulations  and  empirical  measurements  are  described  and  discussed  in  this 
work. 

1.3  Thesis  Outline 

This  thesis  presents  some  background  information  related  to,  the  objectives  of,  and 
the  methodology  to  be  applied  to  my  thesis  research.  It  is  divided  into  the  following 
chapters:  Chapter  II  -  theory  and  background  associated  with  the  problem;  Chapter 
III  -  experiments  performed  during  this  investigation;  Chapter  IV  -  experimental 
results  and  analysis;  and  Chapter  Chapter  V  -  conclusions  and  recommendations  for 
future  work.  Chapter  II  discusses  the  background  information  relevant  to  scintillation 
detection,  neutron  kinematics,  scintillating  fiber  materials  of  interest,  photomultiplier 
tubes,  and  the  concepts  behind  digital  data  acquisition.  Chapter  III  discusses  the 
experimental  and  computational  objectives  to  include  the  methodology,  tools,  and 
equipment  to  be  used  in  accomplishing  them.  Chapter  IV  presents  the  results  of  the 
experiments  performed  during  this  investigation.  Chapter  V  summarizes  the  research 
performed  and  provides  recommendations  for  future  fiber  bundle  work. 
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II.  Theory 


2.1  Overview  of  Scintillators 

Radiation  detection  using  scintillation  light  produced  in  special  materials  is  one  of 
the  oldest  techniques  used  in  nuclear  engineering  today.  Such  scintillation  materials 
should  be  able  to  convert  the  kinetic  energy  of  incident  particles  into  detectable  light 
with  reasonable  efficiency.  The  conversion  of  kinetic  energy  to  light  yield  should  also 
be  linearly  proportional  to  the  energy  deposited  over  a  wide  range.  The  scintillating 
material  should  have  good  optical  characteristics;  it  should  be  transparent  to  the 
wavelength  of  its  own  emission  and  its  index  of  refraction  should  be  close  to  glass  to 
allow  efficient  coupling  of  the  scintillation  light  to  light  sensor  being  used  (typically 
a  PMT). 

There  are  two  types  of  scintillators:  inorganic  and  organic.  Inorganic  scintillators 
in  the  form  of  alkali  halide  crystals,  notably  sodium  iodide  (Nal)  crystals,  are  the 
most  widely  used.  Inorganic  scintillators  typically  have  the  largest  light  output  and 
superior  linearity  but  have  the  disadvantage  of  a  slow  response  time.  The  constituent 
atoms  of  inorganic  materials  typically  have  a  high  density  and  Z- value  which  makes 
them  the  ideal  choice  for  gamma-ray  spectroscopy.  Organic  scintillators  generally  have 
a  faster  response  time  but  a  lower  light  yield  than  inorganic  scintillators.  Organic 
scintillators  are  commonly  used  for  neutron  detection  and  beta  spectroscopy  because 
of  the  prevalence  of  hydrogen  in  organic  materials. 

The  light  produced  during  a  scintillation  event  comes  from  three  mechanisms: 
fluorescence,  phosphorescence,  and  delayed  florescence.  Fluorescence  is  the  prompt 
emission  of  visible  radiation  following  excitation  with  a  decay  approximately  equal  to 
the  lifetime  of  the  an  excited  singlet  state  (between  10-9  and  ICR7  sec).  Phosphores¬ 
cence  involves  a  transition  to  the  excited  triplet  state  with  no  return  to  the  excited 
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singlet  state.  The  decay  in  phosphorescence  is  approximately  equal  to  the  lifetime 
of  the  triplet  state  (between  1CT4  and  10  sec).  Delayed  fluorescence  has  the  same 
emission  spectrum  as  prompt  fluorescence,  but  has  an  emission  time  that  is  only  a 
bit  shorter  than  phosphorescence.  The  delay  is  caused  by  a  brief  transition  from  the 
excited  singlet  state  to  the  triplet  state  and  back  to  the  excited  singlet  state.  A  stan¬ 
dard  way  to  represent  these  quantum  mechanical  state  transitions  is  in  a  Jablonksi 
plot,  as  seen  in  Figure  1.  Scintillators  that  can  convert  a  large  proportion  of  the 


Singlet  Triplet 


Figure  1.  Jablonksi  diagram  of  an  organic  molecule  with  7r— electron  structure  [3]. 

incident  radiation  energy  to  prompt  fluorescence  while  also  minimizing  contributions 
from  phosphorescence  and  delayed  fluorescence  are  typically  used  in  radiation  de- 
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tection.  Scintillators  used  in  pulse  mode  operation  are  primarily  limited  to  prompt 
florescence  because  the  pulse  shaping  times  (for  example,  from  fast  PMTs)  are  much 
shorter  than  typical  phosphorescence  and  delayed  fluorescence  times.  The  photons 
produced  by  these  long-lived  process  are  spread  out  between  signal  pulses  caused  by 
prompt  fluorescence  and  can  be  thought  of  as  noise. 

2.2  Scintillation  Mechanisms  in  Organics 

Organic  scintillator  fluorescence  occurs  from  energy  level  transitions  in  an  organic 
molecule.  For  a  given  molecular  species,  fluorescence  can  be  observed  independent 
of  its  physical  state;  e.g.  a  solid,  a  vapor,  or  part  of  a  solution.  As  a  comparison, 
inorganic  scintillators  like  Nal  require  a  regular  crystalline  lattice  as  a  foundation  for 
the  scintillation  process.  Many  organic  scintillators  are  based  on  organic  molecules 
with  a  7T— electron  structure.  As  seen  in  Figure  1,  energy  can  be  absorbed  by  exciting 
the  electron  configuration  into  any  one  of  the  excited  states.  Singlet  states  (spin  0)  are 
labeled  with  S  and  triplet  (spin  1)  states  are  labeled  with  T.  Molecules  that  are  part 
of  organic  scintillators  have  an  energy  spacing  of  between  3  and  4  eV  from  the  Sq  to 
the  Si  states.  Energy  spacing  between  states  above  S\  being  somewhat  smaller  than 
3  eV.  Each  of  these  electronic  states  are  subdivided  into  energy  states  with  much  finer 
spacing  that  match  various  vibration  states  of  the  molecule.  Spacing  between  these 
vibrational  levels  is  of  the  order  of  0.15  eV,  much  higher  than  average  thermal  energies 
(0.025  eV).  Because  of  this,  practically  all  molecules  are  in  the  lowest  vibrational  state, 
S'oo-  When  an  incident  particle  passing  by  the  organic  scintillator  molecule  deposits 
kinetic  energy,  the  higher  singlet  electronic  states  are  excited  quickly  (on  the  order 
of  picoseconds)  and  then  de-excited  to  the  Si  electron  state  via  internal  conversion. 
Internal  conversion  processes  don’t  emit  any  radiation.  Any  state  that  has  excess 
vibrational  energy  quickly  loses  that  energy  when  reaching  thermal  equilibrium.  The 
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net  effect  of  excitation  is,  after  a  relatively  short  amount  of  time,  a  population  of 
excited  molecules  in  the  S'io  level.  Prompt  fluorescence  scintillation  photons  are 
emitted  in  transitions  between  the  S'io  state  and  the  vibrational  states  of  the  ground 
electronic  state.  If  r  represents  the  fluorescence  decay  time  for  the  S'io  level  (usually 
a  few  nanoseconds  for  organic  scintillators),  the  prompt  fluorescence  intensity  at  time 
t  following  excitation  follows  Equation  1. 

I  =  he~tlr  (1) 

From  Figure  1  it  can  be  seen  how  organic  scintillators  are  transparent  to  their  own 
fluorescence  emission.  Because  all  fluorescence  transitions  have  a  lower  energy  than 
the  minimum  required  for  excitation,  known  as  the  Stokes  shift ,  there  is  little  self- 
absorption  of  the  fluorescence  photons.  Scintillation  efficiency  is  defined  as  the  frac¬ 
tion  of  all  of  incident  particle  energy  converted  into  visible  light.  An  effect  called 
quenching  encompasses  all  de-excitation  processes  that  are  radiation  less.  Impurities 
such  as  dissolved  oxygen  in  liquid  scintillators  can  degrade  the  light  output  by  serving 
as  a  quenching  mechanism  for  excitation  energy. 

2.3  Light  Output  of  Organic  Scintillators 

Most  of  the  kinetic  energy  lost  by  a  particle  incident  upon  a  scintillation  material 
is  dissipated  nonradiatively  in  the  form  of  heat.  The  fraction  of  particle  energy 
converted  to  light,  the  scintillation  efficiency,  depends  on  the  type  of  particle  and 
its  energy.  For  organic  scintillators  the  response  to  electrons  is  linear  up  to  energies 
of  about  125  keV  [4],  Heavy  charged  particles,  such  as  alpha  particles  and  protons, 
always  have  a  smaller  response  when  compared  to  electrons  of  the  same  energy.  The 
response  to  heavy  charged  particles  shows  non-linearity  up  to  larger  energy  thresholds 
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Figure  2.  Energy  distribution  of  recoil  protons  Ep  produced  by  monoenergetic  neutrons. 
For  scatterinf  from  hydrogen,  Ep  is  equivalent  to  Er  [9]. 

when  compared  to  electrons  [9].  For  example,  at  a  few  hundred  keV,  the  response  to 
protons  is  an  order  of  magnitude  smaller  compared  with  the  light  yield  of  equivalent 
energy  electrons  [15].  As  energies  increase  the  discrepancy  is  less,  but  the  proton 
response  is  always  below  the  electron  response.  A  specific  nomenclature  to  used  to 
describe  absolute  light  yield  in  organic  scintillators.  The  MeV  electron  equivalent 
(MeVee)  is  introuced  to  situate  light  yield  on  an  absolute  basis.  The  energy  required 
to  generate  1  MeVee  of  light  by  definition  is  1  MeV  for  fast  electrons  and  several  MeV 
for  heavy  charged  particles  because  of  their  reduced  light  yield  per  energy. 

2.4  Time  Response  of  Organic  Scintillators 

Assuming  that  the  fluorescent  states  in  an  organic  molecule  are  formed  instanta¬ 
neously  and  only  prompt  fluorescence  is  observed,  then  the  time  profile  of  the  light 
pulse  should  be  a  very  fast  leading  edge  followed  by  a  simple  exponential  decay  as 
in  Equation  1.  However,  the  time  it  takes  to  populate  the  fluorescent  states  and  the 
slower  components  of  scintillation  such  as  phosphorescence  and  delayed  fluorescence 
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need  to  be  taken  into  account.  It  takes  approximately  half  a  nanosecond  to  populate 
fluorescent  state  levels  [13].  The  decay  time  from  these  states  is  only  three  to  four 
times  greater  for  very  fast  scintillators  [2],  There  are  two  approaches  that  take  into 
account  this  finite  rise  time.  First,  it  is  assumed  that  the  population  of  the  fluorescent 
states  is  also  exponential  ([10]  and  [14])  and  the  overall  shape  of  the  light  pulse  is 
given  by  Equation  2: 

I  =  /0(e-t/T  —  e-t/Tl)  (2) 

where  r  is  the  time  constant  describing  the  decay  of  the  fluorescent  states  and  T\ 
describes  the  population  of  those  states.  The  second  approach  [16]  to  describing  the 
population  of  the  fluorescent  states  is  to  conclude  that  the  population  step  is  better 
represented  by  a  Gaussian  function  with  a  standard  deviation  get-  The  overall  light 
versus  time  profile  is  then  described  by  Equation  3. 
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The  rise  and  fall  time  of  the  light  output  can  be  experimentally  characterized  by  the 
full  width  at  half  maximum  (FHWM)  of  the  resulting  light  curve.  It  has  become 
commonplace  to  gauge  the  performance  of  very  fast  organic  scintillators  by  their 
FWHM  as  well  as  their  decay  time  [12], 


2.5  BCF-12  Fiber  Scintillators 

Plastic  scintillators  are  organic  scintillators  that  are  dissolved  in  a  solvent  and 
are  polymerized.  A  solvent  consisting  of  a  styrene  monomer  in  which  an  appropriate 
organic  scintillator  is  dissolved  is  a  common  example  of  this  practice  [7].  Other  ex¬ 
amples  of  plastic  matrices  are  polyvinyl  toluene  (PVT)  and  polymethyl  methacrylate 
(PMMA)  [8].  The  fibers  that  will  be  used  for  this  research  are  BCF-12  fibers  from 


Saint-Gobain  Crystals,  Inc.  Their  fibers  consist  of  a  polystyrene  based  core  and  a 
PMMA  cladding.  External  Extra  Mural  Absorber  (EMA)  is  often  used  to  eliminate 
optical  crosstalk.  The  scintillating  core  contains  a  combination  of  fluorescent  dopants 
selected  to  produce  the  desired  scintillation,  optical  and  radiation-resistance  charac¬ 
teristics.  Often,  one  property  is  enhanced  while  another  is  mildly  compromised.  In 
small  fibers  (<0.5mm),  the  fluor  concentration  is  increased,  usually  at  the  expense  of 
light  attenuation  length  [6].  Table  1  lists  the  pertinent  properties  for  the  fibers  used 
in  this  research.  Figure  3  diagrams  what  a  typical  round  scintillating  fiber  is  made 
of. 


Table  1.  BCF-12  Fiber  Properties  [18]. 


Core  Material 

Polystyrene 

Core  Refractive  Index 

1.60 

Density 

1.05 

Cladding  Material 

Acrylic 

Cladding  Refractive  Index 

1.49 

Cladding  Thickness 

3%  of  fiber  diameter 

Numerical  Aperture 

0.58 

Trapping  Efficiency 

3.44%  minimum 

No.  of  H  atoms  per  cc  (core) 

4.82  x  1022 

No.  of  C  atoms  per  cc  (core) 

4.85  x  1022 

No.  of  electrons  per  cc  (core) 

3.4  x  1023 

Radiation  Length 

42  cm 

Emission  Color 

Blue 

Emission  Peak,  nm 

432 

Decay  Time,  ns 

2.7 

1/e  Length,  m 

2.2 

No.  of  Photons  per  MeV 

~  8000 

Operating  Temperature 

-20°  C  to  +50°  C 

Vacuum  Compatible 

Yes 

PMMA  is  the  standard  cladding  material  for  Saint-Gobain  Crystals’  fibers.  It 


has  a  density  of  1.2  g/cc  and  refractive  index  of  1.49. 


The  refractive  index  of  the 
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Figure  3.  Diagram  of  typical  Saint-Gobain  scintillating  fiber  [18]. 
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core  and  cladding  and  the  cross  section  of  the  fiber  determine  the  trapping  efficiency. 
The  trapping  efficiency  also  depends  on  the  distance  between  the  fiber  axis  and  the 
scintillation  event.  For  Saint-Gobain  fibers  the  trapping  efficiency  ranges  from  3.4% 
for  events  occurring  at  the  fiber  axis  to  ~  7%  for  events  near  the  core-cladding 
interface. 

White  or  black  coatings  may  be  applied  to  the  outer  fiber  surface  primarily  to 
eliminate  crosstalk  among  closely  packed  fibers.  Coatings  are  typically  10  to  15 
microns  thick.  An  EMA  coating  decreases  the  overall  signal  intensity  obtained  from 
a  fiber,  irrespective  of  its  length.  This  effect  is  greatest  with  black  EMA,  as  well  as 
with  short  fibers.  The  coating  can  interfere  with  useful  lightpiping  in  the  cladding. 
Black  EMA  applied  at  the  near  end  of  fibers  can  be  used  to  flatten  out  position 
dependent  response.  White  EMA  is  used  in  the  construction  of  short  fiber  imaging 
bundles. 

2.6  Neutron  Kinematics  with  Organic  Scintillators 

One  of  the  most  common  methods  of  fast  neutron  detection  is  based  on  clastic 
scattering  of  neutrons  by  light  nuclei  such  as  hydrogen,  deuterium,  or  helium  [5].  A 
portion  of  the  neutron’s  kinetic  energy  is  transferred  to  the  target  nucleus  resulting 
in  a  recoil  nucleus.  Hydrogen  is  by  far  the  most  popular  target  nuclei  because  it  has 
the  capability  to  transfer  all  of  the  incident  neutron  energy  to  it  as  seen  in  Table  3. 
When  the  target  nuclei  are  light,  the  recoil  nucleus  behaves  much  like  a  proton  or 
alpha  particle  in  the  way  that  it  loses  energy  in  the  detector  medium  [17]. 

The  values  in  Table  3  can  be  derived  from  four  simple  kinematic  equations  [11]. 
These  equations  use  the  variable  definitions  listed  in  Table  2  which  are  graphically 
represented  by  Figure  4.  For  incoming  nonrelativistic  neutrons  (En  <<  939  MeV), 
conserving  momentum  and  energy  in  the  center-of-mass  coordinate  frame  give  the 
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Table  2.  Kinematics  of  neutron  scattering  variables. 


Symbol 

Definition 

A 

mass  of  target  nucleus/neutron  mass 

En 

incoming  neutron  kinetic  energy  in  laboratory  system 

Er 

recoil  nucleus  kinetic  energy  in  laboratory  system 

0 

scattering  angle  of  neutron  in  center-of-mass  coordinate  system 

e 

scattering  angle  of  recoil  nucleus  in  the  laboratory  coordinate  system 

relation  for  the  energy  of  the  recoil  nucleus  given  by 

Er  =  (if^l-cos  @)En.  (4) 

To  convert  to  the  laboratory  coordinate  system,  where  the  target  nucleus  is  at  rest, 
uses  the  transformation 

cos  6  =  sjl- (5) 

Combining  Equations  4  and  5  gives  the  relation  for  the  recoil  nucleus  energy  in  terms 
of  its  own  angle  of  recoil, 

Er  =  (^(cos  29)En.  (6) 

From  Equation  6  it  can  be  seen  that  the  energy  transferred  to  the  recoil  nucleus 
is  determined  by  the  scattering  angle.  If  6  =  90°,  then  the  recoil  energy  is  near  zero. 
However,  a  head-on  collision  between  the  incident  neutron  and  target  nucleus  leads  to 
a  recoil  in  the  same  direction,  0  =  0°,  resulting  in  a  maximum  possible  recoil  energy 
as  seen  in  Equation  7. 


Er 


max 


4A  p 
(1+A)2  n 


(7) 


Table  3  lists  the  maximum  fraction  of  incoming  neutron  energy  that  can  be  trans- 
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ferred  to  a  recoil  nucleus  for  various  target  nuclei.  As  expected,  the  maximum  frac¬ 
tional  energy  transfer  decreases  as  the  nucleus  mass  increases.  This  trend  explains 
why  hydrogen  assumes  the  predominant  role  in  neutron  detection  using  elastic  scat¬ 
tering. 


Table  3.  Maximum  fractional  energy  transfer  in  neutron  elastic  scattering. 


Target  Nucleus 

A 

&R  1 

En  1 maX 

\H 

1 

1 

\H 

2 

0.889 

\He 

3 

0.750 

\He 

4 

0.640 

12  r< 

6  L' 

12 

0.284 

?0 

16 

0.221 

Center— of— mass  system 


Lab  system 


Scattered 


Figure  4.  Neutron  kinematics  [9]. 


The  energy  distribution  of  recoil  nuclei  are  distributed  between  a  minimum  of  zero 
and  a  maximum  value  given  in  Table  3.  Supposing  all  scattering  angles  are  allowed, 
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a  continuum  of  recoil  energies  between  these  two  extremes  is  expected.  If  cr(0)  is 
defined  as  the  differential  cross  section  in  the  center-of-mass  system,  the  probability 
that  the  neutrons  will  be  scattered  into  d@  about  0  is  given  by  Equation  8, 

P(0)d0  =  2vrsin0d0^  (8) 

where  as  is  the  total  scattering  cross  section  integrated  over  all  angles.  However, 
the  distribution  in  recoil  nucleus  energy  is  more  informative  than  the  distribution 
in  recoil  nucleus  scattering  angle.  Because  P(ER)dER  represents  the  probability  of 
creating  a  recoil  with  energy  in  dER  about  ER.  This  leads  to  Equation  9,  which  after 
substituting  into  Equation  8,  gives  Equation  10. 

P(ER)dER  =  P(0)d0  (9) 

P(Er)  =  (10) 

Equation  10  shows  that  the  expected  recoil  energy  continuum  has  the  same  shape 
as  the  differential  cross  section  <x(@)  as  a  function  of  the  center-of-mass  scattering  an¬ 
gle  for  the  incident  neutron.  The  shape  of  <x(@)  tends  to  favor  forward  and  backward 
scattering  for  most  target  nuclei.  For  a  hydrogen  nucleus,  the  scattering  process  is 
isotropic  in  the  center-of-mass  coordinate  system.  The  expected  proton  recoil  energy 
distribution  is  a  rectangle  extending  from  zero  to  the  total  energy  of  the  incident  neu¬ 
tron.  A  detector  based  on  neutron-hydrogen  scattering  should  have  response  function 
that  matches  this  rectangular  shape,  as  seen  in  Figure  2. 
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2.7  Photomultiplier  Tubes 


A  photomultiplier  tube  serves  two  purposes:  the  conversion  of  photons  into  pho¬ 
toelectrons  and  the  multiplication  of  said  photoelectrons.  Conversion  of  photons  into 
photoelectrons  occurs  at  the  photocathode.  The  spectral  response  of  a  PMT  is  de¬ 
termined  by  both  the  PMT  window  material  and  the  photocathode  material.  The 
photocathode  in  a  PMT  converts  incident  light  energy  into  photoelectrons  with  the 
efficiency  of  this  conversion  depending  on  the  wavelength  of  the  incident  light.  The 
sensitivity  of  the  PMT  can  be  documented  using  either  the  quantum  efficiency  or  the 
cathode  radiant  sensitivity.  The  sensitivity  at  long  wavelengths  is  primarily  limited 
by  a  reduction  of  light  absorption  by  the  photocathode  and  the  relatively  low  energy 
imparted  to  the  photoelectron.  At  shorter  wavelengths  the  window  material  becomes 
the  dominant  factor  in  PMT  sensitivity. 

The  time  characteristics  of  a  PMT  are  determined  exclusively  by  the  electron 
trajectories  because  the  time  required  for  photoemission  is  very  short.  The  electron 
transit  time  of  a  PMT  is  the  average  time  difference  between  the  photon  arriving  at 
the  photocathode  and  the  burst  of  electrons  at  the  anode.  An  important  quantity 
when  characterizing  PMT  timing  characteristics  is  the  spread  in  the  electron  transit 
time.  This  value  is  important  because  it  determines  the  time  width  of  the  pulse  of 
electrons  arriving  at  the  anode.  The  distance  between  the  photocathode  and  first 
dynode  is  a  major  factor  in  the  observed  spread  in  transit  time. 

The  shape  of  the  voltage  pulse  produced  at  the  anode  of  a  PM  following  a  scin¬ 
tillation  event  depends  on  the  time  constant  of  the  anode  circuit.  Two  extremes  can 
occur;  one  where  the  time  constant  is  very  large  compared  to  the  decay  time  of  the 
scintillator  and  one  where  the  anode  time  constant  is  much  smaller  than  the  decay 
time  of  the  scintillator.  This  is  illustrated  in  Figure  5. 

When  the  anode  time  constant  is  very  large  you  get  response  from  the  PMT  as 
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Figure  5.  Plots  of  anode  pulse  for  the  two  extremes  of  large  and  small  anode  time 
constant  [9]. 
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seen  at  the  bottom  of  Figure  6.  The  PMT  fully  integrates  all  the  light  from  the 
scintillator  which  comes  out  in  discrete  pulses  as  see  in  the  middle  of  Figure  6.  If  the 
system  response  of  the  PMT  is  well-known,  then  deconvolution  techniques  may  be 
used  to  extract  timing  information  at  the  level  of  individual  scintillation  events. 


time 


I(t) 


-A. 


Figure  6.  Top:  Individual  electrons  leaving  photocathode.  Middle:  Current  at  anode 
after  transit  through  multiplying  structure.  Bottom:  Leading  edge  of  voltage  pulse 
across  anode  circuit  with  long  time  constant.  [9]. 


2.8  Digital  Data  Acquisition 

Digital  data  acquisition  is  a  powerful  technique  that  replaces  many  of  the  pulse- 
shaping  electronics  commonly  used  in  a  nuclear  engineering  laboratory  such  as  a 
preamplifier,  linear  amplifier,  single  channel  analyzer  (SCA),  and  multichannel  an¬ 
alyzer  (MCA).  Raw  pulses  measured  directly  from  the  detector  on  an  oscilloscope 
can  be  stored  digitally  for  further  post-processing.  Once  the  raw  data  is  saved,  the 
behaviors  of  analog  pulse-shaping  electronics  can  be  implemented  via  data  processing 
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algorithms  written  in  computer  programming  languages  such  as  Fortran  or  C++,  as 
well  as  scripting  languages  like  MATLAB  or  Mathematica. 

Digital  data  acquisition  can  bypass  complications  caused  by  connecting  through 
several  different  analog  pulse-shaping  electronics.  For  example,  some  timing  informa¬ 
tion  is  lost  when  measuring  very  fast  events  from  a  detector  as  it  passes  through  a 
preamplifier  and  a  linear  amplifier.  As  the  signal  propagates  through  each  of  these 
pulse-shaping  electronics,  more  and  more  noise  is  added  at  each  step.  This  noise  adds 
to  uncertainty  in  time  measurements.  The  shortest,  least  complicated  path  for  signal 
propagation  will  allow  for  the  most  unambiguous  measurements  from  any  detector. 

In  the  case  of  measuring  scintillation  events,  one  can  measure  directly  from  the 
dynode  or  anode  of  the  PMT  being  used.  This  output  can  then  be  connected  di¬ 
rectly  to  a  fast  digital  oscilloscope  such  as  the  Textronix  DPO7104  as  seen  in  Figure 
7.  Specifications  for  the  Tektronix  DPO7104  can  be  found  in  Tables  4  and  5.  The 
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Figure  7.  Diagram  of  dual-ended  PMT  readout  using  digital  data  acquisition.  The 
light  output  from  the  scintillator  is  converted  to  voltage  pulses  by  the  PMTs,  which 
are  then  measured  by  the  digital  oscilloscope.  The  oscilloscope  captures  a  set  number 
of  waveforms  digitally  to  its  buffer  followed  by  transferring  the  buffer  to  the  computer 
via  a  gigabit  Ethernet  connection. 


18 


DPO7104  has  three  advantages  that  make  it  ideal  for  digitally  acquiring  fast  signal 
pulses.  The  first  advantage  is  its  speed.  It  has  a  maximum  sample  rate  of  20.0  GHz, 
meaning  it  has  a  time  resolution  that  can  get  as  low  at  200  ps  per  datapoint.  The 
second  advantage  is  its  FastFrame  acquisition  capability.  With  FastFrame  acquisi¬ 
tion,  thousands  of  individual  waveforms  per  second  can  be  measured  and  stored  in 
the  oscilloscope’s  memory.  The  only  limitation  is  the  size  of  the  oscilloscope  memory 
and  its  sampling  rate.  The  third  advantage  when  using  a  DPO7104  oscilloscope  is 
the  ability  to  manipulate  every  setting  on  the  oscilloscope  through  Virtual  Instru¬ 
ment  Software  Architecture  (VISA)  commands.  Popular  laboratory  data  acquisition 
programs  like  MATLAB  and  LabView  both  have  built  in  VISA  interfaces  compatible 
with  the  DPO7104.  MATLAB  also  has  instrument  drivers  that  can  be  used  to  directly 
control  instruments  like  the  DPO7104  oscilloscope.  However,  MATLAB  instrument 
drivers  are  not  available  for  every  instrument;  user-created  drivers  can  be  notoriously 
unstable  and  prone  to  data  corruption.  The  VISA  interface  is  independent  of  in¬ 
strument  drivers  and  works  as  long  as  both  the  controller  (e.g.  MATLAB)  and  the 
instrument  (e.g.  Tektronix  DPO7104)  are  VISA  compliant.  An  advantage  to  using 
the  VISA  interface  to  digitally  acquire  data  from  an  oscilloscope,  like  the  DPO7104, 
is  the  capability  to  control  every  feature  of  the  oscilloscope  as  if  one  were  manually 
measuring  signals  on  the  oscilloscope.  For  example,  it  typically  takes  manual  adjust¬ 
ment  of  settings  like  the  voltage  scale,  time  scale,  resolution,  trigger  settings,  and 
waveform  positioning  to  get  the  oscilloscope  to  display  all  of  the  pertinent  features  of 
a  waveform  of  interest.  Once  these  settings  are  “dialed-in”  manually,  the  equivalent 
VISA  commands  can  be  programmed  into  a  MATLAB  script  to  call  up  these  settings 
any  time  they  may  be  needed.  Small  adjustments  like  voltage  scale  or  trigger  level 
can  easily  be  adjusted  this  way  depending  on  changing  experimental  conditions  such 
as  using  a  different  source  or  detector  material. 
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Table  4.  Vertical  acquisition  specifications  for  the  Tektronix  DPO7104  digtal  oscillo¬ 
scope  used  during  this  investigation. 


Input  Channels 

4 

Bandwidth 

1  GHz 

Rise  Time  10%  to  90%  (Typical) 

300  ps 

Rise  Time  20%  to  80%  (Typical) 

200  ps 

DC  Gain  Accuracy 

±1%  with  offset /position  set  to  0 

Bandwidth  Limit 

2.0  GHz 

Input  Coupling 

AC,  DC,  GND 

Input  Impedance 

1  Mil  ±1%  with  13  pF  ±2  pF,  5011 
±1% 

Input  Sensitivity 

1  MG  :  1  mV/div  to  10  V/div 

50  G:  1  mV/div  to  1  V/div 

Vertical  Resolution 

8  bit  (>  11  bit  with  Hi  Res) 

Max  Input  Voltage,  1  MG 

150  V  CAT  I,  derate  at  20 
dB/decade  to  9  Vrms  above  200 
kHz 

Max  Input  Voltage,  5011 

5 Vrms,  with  peaks  <  ±24  V 

Position  Range 

±5  divisions 

Delay  between  any  Two  Channels 

<  100  ps  (50  G) 

The  Tektronix  DPO7104  used  in  this  research  saves  all  waveform  data  internally 
as  signed  8-bit  integers  ranging  from  -128  to  +127.  Then,  depending  on  the  voltage 
scale  chosen,  a  conversion  factor  is  needed  to  convert  from  these  integer  data  values  to 
actual  voltages.  This  highlights  a  disadvantage  of  using  a  digital  oscilloscope  such  as 
the  DPO7104.  The  fidelity  of  voltage  measurements  is  limited  by  the  chosen  voltage 
scale  as  well  as  the  requirement  that  any  pertinent  waveform  data  fit  within  the 
voltage  window  and  time  window  of  the  oscilloscope.  When  taking  measurements 
one  should  keep  in  mind  that  voltage  data  points  will  always  be  separated  by  at  least 
the  voltage  window  range  divided  by  256  (for  8-bit  oscilloscopes). 

To  maximize  the  digital  measurement  capabilities  of  the  oscilloscope,  two  groups 
of  settings  need  to  be  adjusted.  First,  the  baseline  of  a  negative  polarity  pulse  needs 
to  be  set  as  close  to  the  top  of  the  voltage  window  as  possible  without  cutting  off  any 
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Table  5.  Horizontal  acquisition  specifications  for  the  Tektronix  DPO7104  digtal  oscil¬ 
loscope  used  during  this  investigation. 


Maximum  Sample  Rate  (1  ch) 

20GS/s 

Maximum  Sample  Rate  (2  ch) 

lOGS/s 

Maximum  Sample  Rate  (3-4  ch) 

5GS/s 

Maximum  Equivalent  Time  Sam¬ 
pling  Rate 

4  TS/s 

Maximum  Record  Length  with 
Standard  Configuration 

50M  (1  ch),  25M  (2  ch),  12.5M  (3-4 
ch) 

Time  Base  Range 

50  ps/div  to  1000  s/div 

Time  Resolution  (in  ET /IT  mode) 

500  fs 

Time  Base  Delay  Time  Range 

5  ns  to  250  s 

Trigger  Jitter  (RMS) 

1.5  ps rms  with  enhanced  triggering 
OFF 

<100  fs  ft  ms'  with  enhanced  trigger¬ 
ing  ON 

Maximum  FastAcq  Waveform  Cap¬ 
ture  Rate 

>250,000  wfms/s  on  all  4  channels 
simultaneously 

Fast  Frame  Acquisition 

Acquisition  memory  divided  into 
segments;  maximum  trigger  rate 
>310,000  waveforms  per  second. 
Time  of  arrival  recorded  with  each 
event.  Frame  Finder  tool  helps  to 
visually  identify  transients 
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relevant  positive  polarity  portions  of  the  waveform.  This  allows  for  the  entire  range 
of  the  voltage  window  to  be  used  in  measuring  the  pulse  voltage  (using  all  8-bits). 
Second,  a  balance  between  the  time  resolution  and  number  of  data  points  needs  to 
be  struck.  In  other  words,  one  should  keep  the  number  of  data  points  large  enough 
to  describe  the  waveform  being  studied  with  good  statistics  but  not  so  large  as  to  not 
be  able  to  resolve  very  fast  (e.g.  sub-nanosecond)  features  in  the  waveform. 

When  connecting  detectors  directly  to  the  DPO7104  oscilloscope,  an  RG-58  C/U 
BNC  cable  should  be  used  because  of  its  characteristic  impedance  of  50  O.  The 
oscilloscope  can  be  set  to  accept  an  input  impedance  of  1  MO  or  50  O.  The  1  MO 
setting  should  be  used  only  if  the  cable  being  used  is  terminated  with  a  resistor  that 
matches  its  impedance.  For  example,  an  RG-62  cable  would  need  to  be  terminated 
with  a  930  resistor.  Using  50  O  cables  and  setting  the  oscilloscope  input  impedance 
to  50  O  negates  the  need  for  any  terminator  resistors. 
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III.  Methodology 


3.1  Simulation  of  Neutron  Radiation  on  Fiber  Bundles 

A  Monte  Carlo  N-Particle  (MCNP)  transport  model  was  developed  with  the  ob¬ 
jective  of  analyzing  the  interaction  of  neutrons  with  fiber  bundles  of  varying  size. 
Fourteen  fiber  bundles  were  modeled  with  dimensions  as  follows:  lxl,  2x2,  3x3,  4x4, 
5x5,  6x6,  7x7,  8x8,  9x9,  10x10,  15x15,  20x20,  40x40,  and  100x100.  The  fibers  were 
modeled  with  a  diameter  of  500  microns  and  a  fiber  spacing  of  5  microns,  to  match  the 
size  of  the  cladding.  The  space  between  the  fibers,  as  well  as  the  space  encompassing 
the  entire  fiber  bundle  and  source  location,  was  modeled  as  a  perfect  vacuum.  The 
fibers  were  modeled  as  a  bare  styrene  (C6H5CH=CH2)  core  without  the  cladding.  A 
2.5  MeV,  isotropic  neutron  source  was  modeled  and  placed  at  different  locations  rel¬ 
ative  to  the  fiber  bundle.  Five  million  source  particles  were  simulated.  The  different 
source  locations  are  shown  in  Figure  8. 
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Figure  8.  Different  arrangement  of  source  and  fiber  bundle  in  the  MCNP  simulations. 
Parts  A,  B,  and  C  show  the  source  increasing  in  distance  (5,  10,  and  20  cm)  from  the 
fiber  bundle  while  staying  centered  longitudinally  and  vertically  with  the  bundle.  Parts 
D,  E,  and  F  show  the  source  increasing  in  distance  (5,  10,  and  20  cm)  from  the  top-left 
corner  of  the  fiber  bundle  while  staying  centered  longitudinally  with  the  bundle.  Parts 
G  and  H  are  identical  to  Part  A,  except  the  source  is  located  at  the  far  end  and  close 
end  of  the  fiber  bundle  respectively. 
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MATLAB  code  was  written  to  create  the  different  MCNP  input  files  with  some 
variables  entered  by  the  the  user,  to  include:  the  number  of  source  particles,  the  source 
particle  energy,  the  source  location,  the  spacing  between  the  fibers,  and  the  number 
of  fibers  along  the  x-  and  y-axis.  The  fibers  are  always  parallel  to  the  z-axis.  At  the 
2x2  fiber  bundle  size  the  MCNP  input  hie  has  34  total  lines  of  code  -  a  reasonable 
amount  of  lines  to  manually  code  by  hand.  However,  the  100x100  fiber  bundle  input 
hie  has  over  34,000  lines  of  code  -  a  coding  task  better  left  to  a  computer.  Although 
MCNP  has  the  native  capability  to  create  repeated  lattice  structures  such  as  a  bundle 
of  identical  hbers,  this  capability  could  not  be  used  for  two  reasons.  First,  each  hber 
needed  to  be  identified  by  a  unique  cell  number  when  using  the  particle  track  data. 
This  is  not  possible  when  using  a  lattice  framework  in  MCNP.  Second,  the  spacing 
and  material  between  the  hbers  is  more  difficult  to  control  when  using  MCNP’s  lattice 
functionality.  This  second  point  is  important  if  one  should  wish  to  model  a  free-air, 
widely  spaced  hber  bundle.  Once  the  repeating  nature  of  the  hber  bundle  geometry 
is  well  understood,  it  is  relatively  straightforward  to  develop  MATLAB  code  that  can 
create  MCNP  input  hies  that  model  different  sized  hber  bundles  while  still  allowing 
for  each  hber  to  be  individually  addressed.  Run  times  for  creating  the  MCNP  input 
hies  lasted  between  1  and  30  seconds,  while  the  simulation  run  times  lasted  between 
1  to  29  minutes  depending  on  the  size  of  the  bundle  being  processed. 

MCNP  was  configured  to  output  particle  track  information  for  all  particles  that 
entered  the  volume  of  any  hber  in  the  bundle.  MCNP  would  output  the  total  history 
of  these  particles  until  they  left  the  bundle  completely.  Fiber  surface  crossings  and 
collisions  with  the  carbon  or  hydrogen  in  the  hber  were  the  two  types  of  events 
recorded.  At  each  of  these  events,  the  particle’s  position,  direction,  and  energy  were 
recorded.  If  it  was  a  surface-crossing  event,  the  surface  it  crossed  was  also  recorded. 
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For  collision  events,  the  cell  where  the  collision  took  place  and  the  type  of  particle 
collision  (hydrogen  or  carbon)  was  recorded  along  with  the  new  energy  of  the  particle. 

MATLAB  code  was  then  developed  to  parse  through  the  MCNP  ASCII  particle 
track  hies.  These  particle  track  hies  ranged  in  size  from  1.4  MB  for  a  lxl  bundle  to 
8400  MB  for  a  100x100  bundle.  This  code  parses  out  a  list  of  collision  events  saving  the 
type  of  collision  (hydrogen  or  carbon),  the  hber  it  occurred  in,  the  energy  deposited 
for  that  event,  and  the  number  of  total  scattering  events  for  that  an  individual  particle 
history. 

3.2  Characterizing  PMT  Dark  Counts  and  Light  Leakage 

All  PMTs  have  an  intrinsic  dark  current  as  discussed  in  Chapter  II.  For  PMTs 
that  aren’t  sealed  with  a  scintillating  material  during  manufacturing,  there  is  also  the 
potential  for  stray  external  photons  to  reach  the  photocathode  of  the  PMT  if  it  isn’t 
in  a  perfectly  dark  environment.  It  is  important  to  identify  and  characterize  both  the 
dark  current  pulses  and  light  leakage  pulses  so  that  they  can  be  distinguished  from 
the  weak  light  pulses  created  by  organic  scintillators. 

In  order  to  measure  only  the  pulses  arising  from  a  PMT’s  dark  current,  an  abso¬ 
lutely  dark  environment  needs  to  be  created  around  the  PMT’s  window.  The  first 
experiment  to  explore  this  involved  using  an  aluminum  can  and  collar  as  shown  in 
Figure  9.  A  lens  cap  was  placed  in  front  of  the  PMT  window  before  being  sealed  up 
with  the  aluminum  collar  and  can.  Additionally,  electrical  tape  was  triple-wrapped 
around  the  seam  between  the  collar  and  can  as  well  as  around  the  tiny  gap  between 
the  collar  and  PMT  tube  to  prevent  any  stray  light  from  reaching  the  photocathode. 
Next,  one  million  waveforms  were  digitally  acquired  using  the  digital  data  acquisition 
setup  described  in  Chapter  II. 

Another  dark  count  experiment  was  designed  in  order  to  test  the  light-tight  ef- 
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Figure  9.  Aluminum  can  and  collar  system.  The  collar  attaches  to  the  PMT  tube 
via  set  screws.  The  aluminum  can  bolts  to  the  collar.  The  can  holds  either  the  Nal 
crystal  or  the  BC-404  bulk  plastic  scintillator.  For  the  dark  count  experiment,  the  can 
was  empty,  no  scintillating  material  was  present.  The  can  was  used  to  cut  down  light 
leakage  onto  the  PMT.  Tape  is  wrapped  around  the  seam  between  the  collar  and  can. 
More  tape  is  wrapped  around  the  gap  between  the  collar  and  the  PMT  tube. 


fectiveness  of  the  large  dark  box  seen  in  Figure  10.  The  large  interior  space  of  this 
dark  box  lends  itself  well  to  further  experiments,  provided  it  is  sufficiently  light-tight, 
because  it  can  stage  larger  equipment  configurations  and  longer  fiber  bundles.  An 
experiment  combining  the  dual-ended  readout  PMT  setup  with  CCD  cameras  for 
imaging  the  fiber  bundle  is  an  example  of  one  could  take  advantage  of  this  large  in¬ 
terior  space.  A  bare  PMT  was  placed  inside  the  large  dark  box  with  the  high  voltage 
bias  supply  cable  and  signal  cable  connecting  through  the  back  wall  of  the  dark  box 
(as  seen  in  Figure  11)  to  maintain  the  interior  as  light-tight  as  possible.  A  million 
waveforms  were  digitally  acquired  using  this  setup. 

A  similar  setup  to  aluminum  collar  and  can  system  was  used  to  test  whether  or 
not  dark  pulses  can  be  distinguished  from  scintillation  pulses  in  “real-time”  as  they 
are  being  digitally  acquired.  A  cylindrical  piece  of  BC-404  scintillating  plastic  (as 
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Figure  10.  Image  of  the  large  dark  box  used  for  three  experiments  during  this  inves¬ 
tigation.  It  is  designed  to  be  water  tight.  Water  tight  doesn’t  mean  light  tight,  as 
discovered  over  the  course  of  this  research. 


Figure  11.  Left:  Interior  BNC  connections  for  two  pairs  of  PMTs.  Right:  Exterior 
BNC  connections  for  the  PMTs  inside  the  dark  box.  These  connections  allowed  for  the 
PMT  bias  and  signal  cables  to  be  connected  while  keeping  the  box  completely  closed 
and  sealed. 
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seen  on  the  right  side  of  Figure  12)  was  optically  coupled  to  the  PMT,  instead  of 
placing  a  lens  cap  in  front  of  the  PMT  window  -  inside  the  aluminum  can. 


Figure  12.  A  picture  of  the  bulk  scintillator  used  during  this  investigation.  The  Nal 
scintillator  is  on  the  left  and  the  BC-404  plastic  scintillator  is  on  the  right.  Electrical 
tape  is  wrapped  around  the  plastic  scintillator  to  prevent  stray  light  from  radially 
entering  the  scintillator  medium.  The  end  of  the  plastic  scintillator  not  facing  the 
PMT  window  is  also  blacked  out  for  the  same  reason. 


Before  applying  the  optical  grease,  the  PMT  window  was  cleaned  with  isopropyl 
alcohol  and  air-dried.  Bicron  6630-00001  optical  grease  was  applied  thinly  across  the 
face  of  the  scintillating  plastic  before  it  was  coupled  to  the  PMT  window.  Another 
one  million  waveforms  were  digitally  acquired  using  this  setup. 

The  time  necessary  to  acquire  the  one  million  waveforms  in  each  of  the  experi¬ 
ments  just  described  was  recorded  along  with  the  waveforms  themselves.  Because  the 
total  number  of  acquisitions  was  the  same  from  experiment  to  experiment,  a  longer 
acquisition  time  means  the  enclosure  is  doing  a  better  job  keeping  stray  photons 
away  from  the  photocathode  of  the  PMT.  Total  acquisition  time  gives  a  quantitative 
comparison  of  how  light-tight  a  specific  enclosure  was. 


3.3  Construction  of  Fiber  Bundle 


The  material  properties  of  the  BCF-12  plastic  scintillator  fibers  used  in  this  re¬ 
search  can  be  found  in  the  Table  1.  The  first  type  of  fiber  bundle  created  was  a 
5x5  array  with  a  2.54mm  air  gap  between  fibers.  To  start,  25  fibers  were  fed  through 
holes  in  two  layers  of  perfboard  and  cemented  to  the  perfboard  with  steel  impregnated 
epoxy.  The  fibers  were  arranged  vertically  during  the  cementing  step  to  prevent  any 
permanent  bending  of  the  fibers  after  the  epoxy  had  cured.  The  ends  of  the  fibers 
were  cleaved  with  a  wire  cutter.  This  type  of  fiber  bundle  can  be  seen  in  Figure  13. 


Figure  13.  Image  of  first  attempt  at  constructing  a  fiber  bundle.  The  25x25  bundle  is 
framed  by  standard  electronics  perfboard.  The  fibers  are  cemented  to  the  perfboard 
by  a  steel-impregnated  binary  epoxy.  There  is  an  2.54  mm  air  gap  between  each 
fiber.  Bending  of  the  fibers  is  a  noticeable  drawback  of  this  method  of  fiber  bundle 
construction. 


The  first  fiber  bundle  was  abandoned  for  three  reasons.  First,  it  was  difficult 
to  cleanly  apply  the  epoxy  on  the  perfboard  and  between  the  fibers.  The  epoxy  is 
sticky  and  often  ended  up  several  centimeters  along  the  length  of  the  fibers  during 
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the  epoxy  application  step.  Second,  the  method  of  cleaving  the  fiber  ends  left  them 
jagged  -  incapable  of  optically  coupling  to  a  PMT  window.  Third,  even  though  effort 
was  taken  to  keep  the  fibers  taut,  some  of  them  still  managed  to  relax  during  the 
cementing  step  -  leaving  a  slight  bow  over  the  length  of  the  fiber  bundle. 

The  second  type  of  fiber  bundle  created  was  a  7x3  array  with  no  air  gap  between 
the  fibers.  First,  a  ribbon  of  seven  fibers  was  laid  down  over  plastic  supports  spaced 
about  2  inches  apart  from  each  other  as  seen  in  Figure  14.  The  ends  of  the  ribbon  were 


Figure  14.  Image  of  how  the  fiber  ribbons  were  made.  Fibers  are  held  taut  and  clamped 
down  at  both  ends.  They  are  anchored  to  the  four  supports  by  double-sided  tape.  One- 
minute  epoxy  (shown  at  the  top-middle  of  the  image)  is  carefully  applied  to  the  fiber 
ribbon  near  the  clamp  ends  and  in  the  regions  between  the  supports. 

held  taught  by  clamps.  Double-sided  tape  was  used  to  temporarily  fasten  the  ribbon 
to  the  plastic  supports.  After  the  ribbon  was  secured,  Loctite  1-minute  epoxy  was 
applied  with  a  brush  along  the  fiber  bundle  between  the  gaps  of  the  plastic  supports. 
This  epoxy  was  easier  to  work  with  and  completely  cured  in  about  4  hours.  After 
that,  a  second  and  a  third  ribbon  were  made  in  the  exact  same  manner.  Next,  the 
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three  ribbons  were  joined  together.  A  thin  layer  of  epoxy  was  spread  over  the  first 
ribbon  and  the  second  ribbon  placed  on  top.  The  two  ribbons  were  then  clamped 
together  and  left  to  cure  for  about  4  hours.  The  third  ribbon  was  epoxied  to  the  first 
pair  in  the  same  way.  Lastly,  the  fiber  bundle  ends  were  cleaved  with  a  hot  blade. 
The  hot  blade  that  was  used  was  part  of  a  soldering  iron  kit  and  can  be  seen  in  Figure 
15.  The  advantage  of  a  hot  blade  is  its  ability  to  cleave  plastic  fibers  without  the  need 


Figure  15.  Hot  knife  used  to  cleave  ends  of  the  fibers  used  during  this  investigation. 
The  blade  is  adapted  to  attach  to  the  end  of  a  standard  electronics  soldering  iron. 

for  polishing  afterward.  A  hot  blade  simultaneously  polishes  the  fiber  surface  while 
cutting  it  [1].  See  Figures  16  and  17  for  the  end  result  of  the  second  fiber  bundle 
construction  method. 

The  second  fiber  bundle  is  superior  to  the  first  for  three  reasons.  First,  the  fiber 
bundle  is  inherently  more  rigid  because  of  how  the  fibers  are  epoxied  to  each  other. 
Second,  there  is  no  bowing  in  the  second  fiber  bundle  as  there  was  in  the  first.  Third, 
the  fiber  bundle  faces  are  more  reliably  cleaved  and  polished  when  compared  to  the 
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Figure  16.  Face  view  of  the  fiber  bundle  taken  with  CCD  camera  with  magnifying 
optics.  Notice  the  close  packing  of  the  fibers  in  each  row  compared  to  the  epoxy  bulges 
separating  the  three  layers.  Also  note  how  the  top  layer  isn’t  lined  up  with  the  bottom 
two  layers. 


Figure  17.  Image  of  the  fully  constructed  7x3,  no  air-gap  fiber  bundle.  The  bundle  is 
19  cm  long.  A  close-up  of  the  bundle  face  can  be  found  in  Figure  16. 
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first  fiber  bundle  type.  This  will  later  allow  for  optical  coupling  with  the  PMT 
windows. 

In  general  there  are  inherent  complications  when  working  with  plastic  fibers.  One 
complication  arises  from  getting  an  acceptable  cleave  of  the  fiber  end.  Sometimes  the 
hot  blade  would  pinch  more  than  it  would  cut  leading  to  a  gradual  shortening  of  the 
fiber  bundle  length  as  the  perfect  cut  was  sought.  The  individual  fibers  used  in  the 
7x3  bundle  started  out  at  30  cm  while  the  bundle  itself  ended  up  at  19  cm.  Figure  18 
illustrates  the  wasted  bits  of  fiber  left  over  after  finally  getting  an  acceptable  cleave 
of  the  fiber  bundle. 


Figure  18.  Fiber  scraps  from  left  from  botched  attempts  at  cleaving  individual  fibers, 
fiber  ribbons,  and  fiber  bundles.  The  fibers  that  were  used  for  the  fiber  ribbon  con¬ 
struction  started  out  at  30  cm.  When  constructed  the  fiber  ribbons  ended  up  at  25  cm. 
The  final  bundle  measured  in  at  19  cm.  The  difficulty  in  getting  a  clean  fiber  cleave 
resulted  in  wasting  33%  of  the  initial  fiber  length. 
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3.4  Dual-Ended  Readout  of  Fibers 


Measuring  a  scintillation  light  pulse  simultaneously  from  both  ends  of  a  scintil¬ 
lating  fiber  can  potentially  be  used  to  determine  where  the  scintillation  event  took 
place  along  the  length  of  the  fiber.  This  kind  of  dual-ended  readout  is  also  useful  in 
discriminating  scintillation  signals  from  a  noisy  background. 

Before  any  data  collection  using  two  PMTs  can  take  place,  they  need  to  be  gain 
matched.  Gain  matching  ensures  that  two  PMTs  will  produce  the  same  magnitude 
voltage  pulse  from  the  same  amount  of  light  collected.  Six  voltage  pulse  height  spectra 
for  each  PMT  used  in  this  research  were  taken  at  high-voltage  biases  of  2.4,  2.5,  2.6, 
2.7,  2.8,  and  2.9  kV.  The  spectra  were  produced  by  the  scintillation  events  from  a  Nal 
crystal  irradiated  with  a  Cs-137  radioactive  source.  For  each  spectrum,  the  channel 
number  associated  with  the  full  energy  peak  of  the  662  keV  gamma  ray  of  Cs-137 
was  recorded.  As  the  gain  increases,  so  does  the  channel  number  of  the  full  energy 
peak.  Because  a  linear  relationship  exists  between  gain  and  channel  number,  a  least- 
squares  fit  gives  a  linear  equation  for  channel  number  versus  bias  voltage  for  both 
PMTs.  With  these  two  equations,  one  PMT  can  arbitrarily  be  set  to  a  certain  bias, 
while  the  other  is  set  to  the  bias  that  will  produce  the  same  response  seen  in  the  first 
PMT. 

Three  different  dual-ended  readout  techniques  were  attempted  during  this  investi¬ 
gation,  although  only  one  was  used  for  actual  data  collection.  The  first  and  simplest 
technique  used  two  bare  PMTs  connected  inside  the  large  dark  box  seen  in  Figure  10. 
The  fiber  bundle  was  sandwiched  between  the  two  PMT  windows  and  both  ends  of 
the  bundle  were  coated  with  optical  grease  to  enhance  the  optical  coupling  between 
the  fiber  faces  and  PMT  windows.  This  technique  seemed  to  work  well  for  the  7x3 
fiber  bundle  but  not  for  the  lone  fiber  or  the  small  2x1  bundle.  The  7x3  bundle 
had  a  natural  rigidity  that  kept  the  fibers  from  bending  when  pinched  between  the 
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PMTs.  The  single  fiber  and  2x1  bundle  tended  to  slip  along  the  PMT  window  during 
handling  of  the  dark  box.  This  slipping  is  illustrated  in  Figure  19.  It  was  discovered 
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Figure  19.  Illustration  of  fiber  slipping  along  the  PMT  window.  If  the  fibers  and  PMTs 
are  not  anchored  securely,  then  fiber  slipping  can  occur.  Fiber  slipping  is  problematic 
for  two  reasons:  the  fiber  face  moves  away  from  the  center  (most  sensitive  part)  of 
the  PMT  window  and  light  from  the  fiber  face  has  a  chance  to  skip  across  the  PMT 
window  without  reaching  the  photocathode. 


during  this  attempt  that  both  PMTs  needed  to  be  securely  anchored  to  the  dark  box 
(after  just  enough  compression  was  applied  to  the  fiber  bundle)  to  ensure  the  fibers 
did  not  slip  during  measurement. 

This  led  to  the  second  attempt  at  reliably  acquiring  dual-ended  readout  data.  A 
cardboard  box  was  triple-wrapped  in  electrical  tape  with  an  opening  flap  on  top  and 
circular  holes  cut  on  both  ends  of  the  box  where  PMTs  could  poke  through.  First, 
one  PMT  was  secured  in  place  with  electrical  tape  so  that  its  window  was  inside 
the  cardboard  dark  box.  Second,  the  fiber  bundle  was  placed  on  a  plastic  stand 
inside  the  box  and  coupled  with  optical  grease  to  the  PMT.  Next,  the  remaining 
PMT  was  inserted  through  the  other  end  of  the  box,  pinching  the  fiber  bundle,  and 
secured  with  electrical  tape.  The  opening  flap  was  sealed  and  the  entire  apparatus 
was  placed  inside  the  larger  dark  box  as  seen  in  Figure  20.  This  method  temporarily 
solved  the  issue  of  fiber  slipping  until  the  temperature  around  the  cardboard  dark  box 
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Figure  20.  Image  of  the  cardboard  dark  box.  Electrical  tape  is  triple-wrapped  around 
all  faces  and  edges  of  the  box.  Two  holes  are  cut  on  opposite  ends  of  the  box  to  allow 
for  insertion  of  the  PMTs.  A  hatch  is  cut  on  the  top  of  the  box  to  allow  for  fiber  bundle 
and  radioactive  source  placement.  The  advantage  of  this  design  is  it’s  darker  than  the 
larger  plastic  dark  box  used  previously,  while  also  being  portable.  The  disadvantage 
of  this  design  is  the  fragile  nature  of  the  tape  adhesive.  Changes  in  temperature  made 
the  tape  peel  off  and  loosen  the  PMTs  from  the  box,  resulting  in  fiber  slipping. 


increased  beyond  87°F,  the  upper  temperature  limit  of  the  tape’s  adhesive.  This  rise 
in  temperature  was  caused  by  a  combination  of  the  heat  created  by  the  PMTs  and  an 
increase  in  laboratory  temperature  caused  by  a  malfunctioning  building  heating  and 
cooling  system.  The  heat  caused  the  tape  to  peel  back  and  loosen  the  PMTs  from 
the  cardboard  dark  box.  This  caused  the  PMTs  to  move  and  allow  the  fiber  bundle 
face  to  slip  again. 

The  fragile  nature  of  the  cardboard  dark  box  led  to  the  design  of  an  aluminum 
dark  box  with  similar  functionality.  Six  1/4”  thick  aluminum  sheets  were  cut  to  form 
the  sides  of  the  dark  box.  Two  holes  were  cut  to  the  diameter  of  the  PMT  tubes  on 
opposite  sides  of  the  box.  Two  aluminum  collars  were  also  made  to  mount  to  these 
sides  of  the  aluminum  dark  box  and  secure  the  PMTs  with  set  screws.  The  walls  of 
the  dark  box  were  then  bolted  together  and  electrical  tape  was  triple  wrapped  along 
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all  seams  except  for  the  seam  along  the  lid  of  the  box.  The  first  PMT  was  set  in  place 
and  locked  down  with  set  screws.  The  fiber  bundle  was  placed  on  its  stand  inside 
the  aluminum  dark  box  and  optically  coupled  to  the  first  PMT.  The  second  PMT 
was  inserted  into  the  box  so  that  a  small  amount  of  compression  pinched  the  fiber 
bundle  between  the  two  PMTs.  Once  compression  was  established,  the  second  PMT 
was  secured  with  set  screws.  Figure  21  shows  what  the  fully  constructed  aluminum 
dark  box  looks  like.  Reusable  putty  was  placed  along  the  gap  between  the  PMTs 


Figure  21.  Full  view  of  the  aluminum  dark  box  used  for  data  collection  during  the 
dual-ended  readout  portion  of  this  investigation.  The  PMTs  are  secured  to  the  box 
with  set  screw  after  pinching  the  fiber  bundle  inside  the  box.  The  box  is  designed  with 
portability  in  mind,  being  placed  in  front  of  radiation  sources  that  are  not  found  in  the 
laboratory. 


and  the  aluminum  collars  to  prevent  light  from  getting  inside  the  box.  Putty  was 
used  because  it  could  fill  the  gap  and  could  also  be  easily  moved  whenever  the  PMT 
needed  to  be  moved.  See  Figure  22  for  what  the  putty  looks  like  around  the  PMT 
tube.  The  lid  was  bolted  down  and  triple-sealed  with  electrical  tape  after  the  fiber 
sample  was  secured  and  the  radioactive  source  of  interest  was  placed  inside  the  box. 
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Figure  22.  Close-up  view  of  the  PMT  side  of  the  aluminum  dark  box.  Note  the 
electrical  tape  wrapped  around  all  seams  of  the  box.  Also  worth  noting  is  the  blue 
putty  wrapped  around  the  PMT  tube  near  the  aluminum  collar.  The  putty  is  squeezed 
between  the  gap  formed  by  the  aluminum  collar  and  PMT  tube.  It  was  used  because 
it  is  easy  to  pick  up  and  move  whenever  the  PMTs  needed  to  be  moved. 


Three  radioactive  sources  were  used  during  this  investigation:  Cs-137  for  gamma 
radiation,  Arn-241  for  alpha  radiation,  and  Sr-90  for  beta  radiation.  Table  6  lists 
pertinent  source  information.  The  lone  fiber,  2x1  bundle,  and  7x3  bundle  were 


Table  6.  Radioactive  sources  used  during  dual-ended  readout  portions  of  this  research. 


Isotope 

Radiation 

Energy  [MeV] 

Activity  [//  Ci] 

Cs-137 

Arn-241 

Sr-90 

Gamma 

Alpha 

Beta 

0.662 

5.64 

0.546  (max) 

7.48 

0.00019 

0.074 

separately  irradiated  with  all  three  sources.  Each  source  was  placed  at  three  different 
positions  along  the  fiber  bundle:  1  cm  from  the  left  end,  the  center  of  the  bundle, 
and  1  cm  from  the  right  end  of  the  bundle.  All  three  sources  were  uncollimated  and 
placed  directly  on  the  fiber.  An  additional  experiment  was  performed  using  the  Cs- 


38 


137  source.  It  was  placed  along  the  fibers  in  the  same  manner  described  previously, 
but  the  gammas  were  collimated  using  a  lead  block  and  the  source  was  located  4  cm 
below  the  fiber  bundle. 


3.5  Exploration  of  CCD  Capabilities 

An  investigation  into  the  feasibility  of  using  a  commercial  off-the-shelf  (COTS) 
CCD  camera  to  capture  scintillation  events  coming  from  the  fiber  bundle  was  con¬ 
ducted.  It  is  hoped  that  by  combining  timing  information  from  the  PMTs  and  position 
information  from  a  CCD  image  of  lit  fibers,  one  can  determine  the  path  and  energy  of 
incident  radiation  through  the  bundle.  The  specifications  of  the  camera  used  during 
this  investigation  can  be  found  in  Table  7. 


Table  7.  Specifications  of  the  CCD  camera  used  during  this  investigation. 


Resolution 

1024(H)  x  768(V)  Color  or  Mono 

Pixel  Pitch 

4.65  /mi  x  4.65  /im 

Active  Area 

4.76  mm  x  3.57  mm  -  5.95  mm  di¬ 
agonal 

Max  Datarate 

20  MHz 

Responsivity 

16.2  -  63.9  DN/(nJ/cm2) 

Dynamic  Range 

59.7  dB 

Bit  Depth 

8  and  12-bit 

Exposure  Range 

63.5  {is  to  10  seconds 

Gain 

0  dB  to  24.57  dB  in  29  increments 

Two  experiments  were  designed  to  test  the  feasibility  of  imaging  scintillation 
events  with  a  CCD  camera.  Both  experiments  take  place  inside  the  large  dark  box 
seen  in  Figure  10.  The  setup  for  the  first  experiment  uses  a  Nal  crystal  placed  at  the 
focal  length  of  the  CCD  camera.  Half  of  the  Nal  crystal  was  masked  by  electrical 
tape,  leaving  the  other  half  bare.  A  Cs-137  radioactive  source  was  placed  directly  on 
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top  of  the  Nal  crystal.  The  purpose  of  the  mask  is  to  create  two  distinct  regions  in 
the  camera’s  held  of  view  (FOV).  One  region  should  always  be  dark,  while  the  other 
region  can  transmit  the  scintillation  light  to  the  camera.  The  location  of  the  line 
separating  these  two  regions  is  known  beforehand  so  that  comparisons  can  be  made 
when  images  are  taken  of  the  scintillator  in  the  dark.  Figure  23  shows  the  setup  for 
this  experiment. 


Figure  23.  Experimental  setup  of  Nal  crystal  imaging  with  the  CCD  camera  inside  the 
large  dark  box.  The  Nal  crystal  is  placed  at  the  focal  legnth  of  the  camera  optics.  Half 
of  the  crystal  is  covered  with  tape.  The  camera’s  FOV  captures  both  the  mask  and  the 
exposed  crystal  face.  A  Cs-137  source  is  placed  on  top  of  the  crystal.  Any  scintillation 
light  captured  by  the  camera  will  only  show  up  on  the  right  side  of  image. 


The  second  experiment  places  the  front  face  of  the  7x3  fiber  bundle  at  the  focal 
length  of  the  CCD  camera  system.  Again,  the  Cs-137  source  was  placed  directly  on 
top  of  the  bundle.  The  location  of  each  fiber  in  the  bundle  was  imaged  beforehand 
so  that  comparisons  could  be  made  when  the  scintillator  images  are  taken  in  dark. 
Figure  24  shows  the  setup  for  this  experiment. 
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Figure  24.  Experimental  setup  of  7x3  fiber  bundle  imaging  with  the  CCD  camera  inside 
the  large  dark  box.  The  closest  fiber  bundle  face  to  the  camera  is  located  at  the  focal 
length  of  the  camera  optics.  A  Cs-137  source  is  placed  on  top  of  the  bundle,  centered 
along  the  length  of  the  fibers. 
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IV.  Analysis  of  Results 


4.1  Simulation  of  Neutron  Radiation  on  Fiber  Bundles 

The  results  of  the  MCNP  simulations  are  presented  in  this  section.  A  total  of  42 
simulations  were  run,  split  between  14  simulations  for  each  of  the  four  different  source 
locations.  The  14  simulations  simulated  lxl,  2x2,  3x3,  4x4,  5x5,  6x6,  7x7,  8x8,  9x9, 
10x10,  15x15,  20x20,  40x40,  and  100x100  fiber  bundles.  Particle  track  information 
was  saved  for  each  simulation.  The  model  was  set  up  as  described  in  Chapter  III. 
These  files  were  parsed  using  the  MATLAB  script  located  in  Appendix  E.  The  run 
time  for  processing  these  particle  track  files  ranged  from  about  30  seconds  for  the  lxl 
bundle  simulations  to  over  5  hours  for  the  100x100  bundle  simulations. 

The  energy  deposition  spectra  for  all  fiber  bundles  simulated  show  the  same  basic 
features  as  seen  in  Figure  25.  Note  the  full  energy  peak  that  begins  to  appear  at 
2.5  MeV  from  the  15x15  bundle  onward,  matching  the  incident  neutron  energy  of  2.5 
MeV.  The  peak  occurring  at  around  0.70  MeV  corresponds  to  collisions  with  carbon 
in  the  plastic  scintillator.  This  matches  with  the  theoretical  maximum  values  of  frac¬ 
tional  energy  transfer  a  neutron  can  transfer  during  clastic  scattering  as  documented 
in  Table  3  of  Chapter  II. 

The  results  of  the  simulations  where  the  source  was  located  5  cm  from  the  left 
side  of  the  bundle  while  staying  centered  longitudinally  and  vertically  can  be  found 
in  Figure  26.  As  can  be  seen  in  the  figure,  a  gradient  in  energy  deposition  starts  to  be 
seen  starting  at  the  10x10  fiber  bundle  onward.  As  the  fiber  bundle  increases  in  size, 
so  does  the  intrinsic  efficiency  as  seen  in  Figure  27.  A  40x40  bundle  has  6.7  ±  0.2% 
intrinsic  efficiency  for  double  scatter  events. 

The  results  of  the  simulations  where  the  source  was  located  5  cm  from  the  upper 
left  corner  of  the  bundle,  while  staying  centered  both  longitudinally  and  vertically, 
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Figure  25.  Progression  of  energy  deposition  spectra  for  lxl,  2x2,  3x3,  4x4,  5x5,  6x6, 
7x7,  8x8,  9x9,  10x10,  15x15,  20x20,  40x40,  and  100x100  fiber  bundles  (left  to  right, 
top  to  bottom).  The  incident  neutron  energy  is  2.5  MeV.  Note  the  full  energy  peak 
that  begins  to  appear  at  2.5  MeV  from  the  15x15  bundle  onward.  Also  worth  noting 
is  the  peak  occurring  at  around  0.70  MeV  corresponding  to  collisions  with  carbon  in 
the  plastic  scinillator.  These  features  are  common  to  all  of  the  plots  and  become  more 
apparent  with  larger  bundle  size. 
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Figure  26.  MCNP  simulated  energy  deposited  per  fiber  in  a  lxl,  2x2,  3x3,  4x4,  5x5, 
6x6,  7x7,  8x8,  9x9,  10x10,  15x15,  20x20,  40x40,  and  100x100  fiber  bundle  (left  to 
right,  top  to  bottom).  Lighter  colors  represent  more  energy  deposited  than  darker 
colors.  The  location  of  the  isotropic  neutron  source  is  defined  by  Part  A  of  Figure  8. 
A  noticeable  gradient  in  energy  deposition  can  be  seen  starting  from  the  15x15  bundle 
onward. 
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Figure  27.  Intrinsic  scattering  efficiencies  for  single,  double,  and  triple  scatter  events 
as  a  function  of  bundle  size. 


can  be  found  in  Figure  40.  As  seen  previously,  a  gradient  in  energy  deposition  starts 
to  be  seen  starting  at  the  10x10  fiber  bundle  onward.  One  can  clearly  observe  the 
source  position  has  changed  between  this  group  of  simulations  and  the  ones  discussed 
earlier  by  the  direction  of  this  gradient.  This  energy  deposited  is  again  sufficient  to 
produce  measurable  light  output  from  each  fiber. 

The  results  of  the  simulations  where  the  source  was  located  5  cm  from  the  left 
side  of  the  bundle  while  staying  centered  vertically,  but  located  at  the  far  end  of  the 
fiber  bundle  can  be  found  in  Figure  41  in  Appendix  A.  Figure  42  in  Appendix  A  has 
the  same  source  parameters,  except  the  source  is  at  the  near  end  of  the  fiber  bundle 
instead  of  the  far  end. 

An  estimate  of  the  time  separation  between  double-scatter  events  was  determined 
using  the  particle  track  results  from  40x40  bundle  with  the  source  closest  to  the 
bundle.  First,  the  distances  between  each  double-scatter  event  were  calculated.  The 
distribution  of  these  scatter  events  can  be  found  in  Figure  [?].  The  kinetic  energy  for 
each  of  these  double-scatter  neutrons  during  their  transit  from  the  first  scatter  site  to 
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the  second  scatter  site  was  also  calculated.  The  distribution  of  neutron  energies  can 
be  found  in  Figure  [?].  Using  non-relativistic  kinematic  equations,  the  distance  from 
the  first  scatter  site  to  the  second  scatter  site,  and  neutron  energy  during  transit, 
the  time  separation  between  double-scatter  events  was  calculated.  The  distribution 
of  double-scatter  time  separations  can  be  found  in  Figure  [?]. 


40x40  Bundle  Double-Scatter  Distances 


Distance  [cm] 

Figure  28.  Distribution  of  scatter  distances  for  double-scatter  events  in  the  40x40  fiber 
bundle. 


4.2  Characterizing  PMT  Dark  Counts  and  Light  Leakage 

The  results  of  characterizing  the  dark  counts  of  the  PMTs  used  during  this  inves¬ 
tigation  are  presented  in  this  section.  Three  methods  at  creating  a  dark  environment 
around  the  PMTs  were  studied.  These  results  of  these  methods  are  summarized  in 
Table  9.  As  expected,  the  lowest  dark  count  rate  was  achieved  with  the  lens  cap  / 
aluminum  can  and  collar  system  with  a  dark  count  rate  of  15  per  second.  This  setup 
is  unrealistic  for  two  reasons:  there  is  no  scintillator,  hence  no  radiation  detection 
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Figure  29.  Distribution  of  transit  energies  for  double-scatter  events  in  the  40x40  fiber 
bundle. 


Time  [ns] 

Figure  30.  Distribution  of  scatter  times  for  double-scatter  events  in  the  40x40  fiber 
bundle.  Over  40%  of  the  double  scatter  events  occur  with  a  time  separtion  of  1  ns  or 
greater. 


Table  8.  Dark  current  and  light  leakage  counts  for  the  three  dark  configurations  used 


in  this  investigation. 


Time  to  Million  Counts  [h] 

Rate  [Dark  Counts/s] 

Aluminum  Can  /w  Lens  Cap 

18.7 

15 

Bare  PMT  in  Large  Dark  Box 

0.68 

406 

Aluminum  Dark  Box 

1.2 

229 
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can  take  place,  and  a  cap  completely  covered  the  PMT  window.  The  benefit  of  these 
results  is  they  give  a  baseline  of  how  many  dark  counts  should  be  expected  in  an  ab¬ 
solutely  dark  environment.  The  other  two  methods  produced  dark  count  rates  that 
were  26x  (for  the  PMT  in  the  large  dark  box)  and  16x  greater  (in  the  aluminum  dark 
box).  However,  the  setup  of  both  of  these  experiments  is  representative  of  how  data 
needs  to  be  collected  for  dual-ended  readout  experiments. 

Since  a  noisy  environment  is  expected  and  the  scintillation  pulses  may  be  as  weak 
as  dark  count  pulses,  a  method  to  distinguish  between  these  two  types  of  pulses  is 
required.  The  experiment  using  the  BC-404  bulk  plastic  scintillator  addresses  this. 
It  is  assumed  that  dark  counts  measured  from  a  PMT  arise  from  the  dark  current 
intrinsic  to  the  PMT  or  from  a  few  photons  of  stray  light  hitting  the  photocathode. 
As  discussed  in  Chapter  II,  the  light  produced  during  scintillation  is  not  produced 
instantaneously,  but  instead  has  a  time-dependence.  Therefore,  a  scintillation  pulse 
should  have  a  measurably  wider  width  than  a  pulse  caused  by  a  dark  event.  An 
example  of  a  dark  count  waveform  can  be  found  in  Figure  31,  while  an  example  of  a 
scintillation  event  waveform  can  be  found  in  Figure  32. 

Using  the  assumptions  about  voltage  pulses  from  the  output  of  a  PMT  discussed  in 
Chapter  II,  an  algorithm  was  written  in  MATLAB  to  calculate  the  decay  constant  for 
each  waveform.  A  bimodal  distribution  of  decay  constants  was  discovered,  illustrated 
in  Figure  33.  These  two  distributions  are  statistically  different  from  each  other  using 
the  standard  3cr  test  as  seen  from  the  data  presented  in  Table  9.  The  shorter  decay 
constant  distribution  represents  pulses  arising  from  PMT  dark  current  or  light  leakage. 
The  longer  decay  constant  distribution  comes  from  scintillation  events  in  the  bulk 
plastic  scintillator.  Boundary  values  that  delineate  dark  pulses  from  scintillation 
pulses  can  be  determined,  once  the  distribution  of  decay  constants  is  known.  A  check 
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Figure  31.  An  example  of  a  dark  count  waveform  captured  via  a  PMT  made  light-tight 
with  the  aluminum  can  and  collar  system  seen  in  Figure  9. 


Figure  32.  An  example  of  a  bulk  plastic  scintillation  event  waveform  captured  via  a 
PMT  made  light-tight  with  the  aluminum  can  and  collar  system  seen  in  Figure  9. 


Table  9.  Dark  current  and  light  leakage  counts  for  the  three  dark  configurations  used 


in  this  investigation. 


Time  to  Million  Counts  [h] 

Rate  [Dark  Counts/s] 

Aluminum  Can  /w  Lens  Cap 

18.7 

15 

Bare  PMT  in  Large  Dark  Box 

0.68 

406 

Aluminum  Dark  Box 

1.2 

229 
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Figure  33.  Comparison  of  the  distribution  of  decay  constants  for  PMT  dark  counts 
and  bulk  plastic  scintillation  events.  The  distributions  are  separated  from  each  other 
by  3  standard  deviations.  Discrimination  of  noise  events  from  scintillation  events  can 
be  done  at  acquisition  time  using  these  distributions. 


for  these  values  can  be  inserted  directly  into  the  digital  data  acquisition  code  for 
real-time  discrimination  of  dark  events  from  scintillation  events  using  a  single  PMT. 


4.3  Dual-Ended  Readout  of  Fibers 

The  results  of  the  dual-ended  readout  experiments  are  presented  in  this  section. 
As  discussed  in  Chapter  111,  the  two  PMTs  used  in  this  investigation  needed  to  be  gain 
matched  before  any  data  was  collected.  Gain  matching  ensures  that  two  PMTs  will 
produce  the  same  magnitude  voltage  pulse  from  the  same  amount  of  light  collected. 
Slight  differences  in  the  materials  and  construction  can  cause  a  different  response  to 
the  same  high  voltage  bias  for  two  different  PMTs.  The  data  plotted  in  Figure  34 
uses  the  gain  matching  procedure  outlined  in  Chapter  III.  By  setting  the  two  lines  in 
those  plots  equal  to  each  other,  it  was  determined  that  if  PMT  A  was  set  to  2.70  kV, 
then  PMT  B  would  need  to  be  set  to  2.74  kV  to  provide  the  same  response  as  PMT 
A. 
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Summary  of  Fit 

RSquare  0  995357 

RSquare  Adj  0  994196 

Root  Mean  square  Error  0.12 

Mean  o(  Response  5  1 

Observations  (or  Sum  Wgts)  6 

Analysis  of  Variance 


Sum  of 

Source  OF  Squares  Mean  Square  F  Ratio 

Model  1  12  348000  12  3480  857  5000 

Error  4  0  057600  0.0144  Proto  >F 

C  Total  5  12  405600  <0001' 

Parameter  Estimates 

Term  Estimate  Std  Error  t  Ratio  Proto>|t| 

intercept  -1716  0.761742  -22.53  <0001* 

Bias  Voltage  8  4  0  286855  2  9  28  <  0001* 


Bivariate  Fit  of  PMT  B  By  Bias  Voltage 
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Sum  of 

Source  OF  Squares  Mean  Square  F  Ratio 

Model  1  10  522566  10  5226  682  9034 

Error  4  0  061634  0.0154  Prob>F 

C  Total  5  10  584200  <  0001- 

Parameter  Estimates 

Term  Estimate  Std  Error  t  Ratio  Prob>|i| 

intercept  -15  74386  0  787967  -19  99  <0001* 

Bias  Voltage  7  7542857  0  29673  26  13  <0001* 


Figure  34.  Linear  regression  analysis  used  to  gain  match  the  two  PMTs  used  during 
this  investigation.  From  this  data  it  was  determined  that  PMT  A  would  be  set  to  a  bias 
of  2.70  kV,  while  PMT  would  be  set  to  a  bias  of  2.74  kV.  These  bias  settings  ensured 
the  response  of  either  PMT  would  be  the  same  for  an  identical  scintillation  light  pulse. 
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Data  from  dual-ended  readouts  can  be  analyzed  to  discriminate  out  scintillator 
events  from  background  noise.  The  raw  mean  voltage  distributions  without  noise 
discrimination  can  be  seen  in  Figure  35.  It  can  be  seen  from  that  distribution  that 


Figure  35.  Mean  voltage  distribution  for  an  Am-241  source  placed  at  the  center  of  a 
single  fiber.  The  distribution  is  actually  comprised  to  two  separate  distributions.  The 
peak  at  the  far  left  shows  contribution  from  noise  events  such  as  dark  current  pulses 
or  light  leakage  events.  The  broader  distribution  to  the  right  of  that  peak  represents 
actual  scintillation  events. 


the  noise  events  would  skew  the  mean  value  of  the  scintillation  pulses  toward  lower 
voltages.  Scintillation  events  are  measured  as  pulses  that  occur  within  the  same 
time  window  on  both  PMTs.  Typically,  non-scintillation  events  like  dark  current 
pulses  or  light  leakage  photons  will  only  be  measured  on  one  PMT,  but  not  on  the 
other  PMT  within  the  same  time  window.  An  algorithm  was  developed  to  determine 
the  distribution  of  mean  voltages  measured  from  PMT  A  whenever  no  pulse  from 
PMT  B  rose  3  standard  deviations  above  the  noise  floor.  This  algorithm  determined 
that  the  noise  events  had  a  maximum  mean  voltage  of  0.6  mV.  This  value  was  then 
discriminated  out  of  the  distributions  of  mean  voltages  for  the  single  fiber  and  7x3 
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fiber  bundle  irradiated  by  Cs-137,  Am-241,  and  Sr-90.  Table  10  lists  the  number  of 
scintillation  events,  the  mean  voltage  read  from  PMT  A  and  the  mean  voltage  read 
from  PMT  B  for  the  single  fiber  irradiated  by  these  sources.  Table  11  lists  the  same 


Table  10.  Mean  voltage  pulse  values  for  scintillation  events  in  a  single  BCF-12  scintil¬ 
lating  fiber. 


Isotope 

Scintillation  Events 

Ch  1  Mean  Voltage  [mV] 

Ch  2  Mean  Voltage  [mV] 

Cs-137 

8,477 

2.2 

1.9 

Am-241 

29,697 

1.4 

1.2 

Sr-90 

15,520 

1.4 

1.2 

information,  except  for  the  7x3  bundle  instead  of  the  single  fiber. 


Table  11.  Mean  voltage  pulse  values  for  scintillation  events  in  a  7x3  bundle  of  BCF-12 
scintillating  fibers.  Note  the  increase  in  mean  voltages  compared  to  the  values  in  Table 
10 


Isotope 

Scintillation  Events 

Ch  1  Mean  Voltage  [mV] 

Ch  2  Mean  Voltage  [mV] 

Cs-137 

271,181 

2.2 

1.9 

Am-241 

210,038 

1.5 

1.2 

Sr-90 

433,243 

2.2 

1.9 

As  expected,  pulses  measured  from  the  7x3  bundle  are  larger  than  those  measured 
from  the  single  fiber.  This  is  because  the  interaction  medium  is  larger  in  the  fiber 
bundle;  the  light  produced  by  multiple  fibers  should  be  greater  than  that  produced 
by  a  single  fiber.  Discriminating  out  the  noise  events  allows  only  the  scintillation 
waveforms  to  be  analyzed.  A  representative  waveform  measured  from  the  single  fiber 
irradiated  with  Sr-90  is  seen  in  Figure  37.  Compare  this  with  to  the  7x3  bundle  irra¬ 
diated  with  the  same  source,  as  seen  in  Figure  ??.  Notice  the  fiber  bundle  waveform 
is  larger  than  the  single  fiber  waveform.  There  are  also  more  peaks  seen  in  the  fiber 
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Figure  36.  Representative  waveform  measured  from  the  single  fiber  irradiated  with 
Sr-90. 


bundle  waveforms  compared  to  the  single  fiber.  This  could  indicate  that  the  PMT  is 
fast  enough  to  see  individual  fibers  lighting  up  in  the  bundle.  Representative  wave¬ 
forms  for  single  fibers  and  fiber  bundles  irradiated  with  Am-241  and  Cs-137  can  be 
found  in  Appendix  B. 

The  expected  scintillation  pulses  measured  from  a  PMT  are  characterized  by  a 
fast  rise  time  followed  by  a  slower  exponential  decay,  ffowever,  the  pulses  measured 
from  the  single  fibers  and  the  fiber  bundles  show  multiple  peaks  instead  of  a  single 
decaying  exponential.  The  timing  separation  between  these  peaks  as  well  as  the  total 
number  of  peaks  were  studied.  For  the  purposes  of  this  analysis,  a  peak  was  defined 
as  a  local  maxima  that  was  at  least  6cr  above  the  noise.  The  MATLAB  code  used  to 
locate  these  peaks  can  be  found  in  Appendix  G.  The  distributions  of  time  separation 
between  the  peaks  of  individual  scintillation  pulses  for  single  fiber  and  a  fiber  bundle 
irradiated  with  Sr-90  can  be  found  in  Figure  38.  The  distributions  of  the  number 
of  peaks  measured  is  found  in  Figure  39  If  the  multiple  peaks  were  characteristic 
of  a  single  fiber  firing,  then  multiple  scatter  events  in  the  fiber  bundle  should  have 
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Figure  37.  Representative  waveform  measured  from  the  7x3  fiber  bundle  irradiated 
with  Sr-90. 


Figure  38.  Distributions  of  peak  separation  times  measured  during  Sr-90  irradiation  of 
single  fiber  and  7x3  fiber  bundle. 


55 


10 


Sr-90  Scintillation  Peak  Distribution 
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Figure  39.  Distributions  of  number  of  peaks  measured  during  Sr-90  irradiation  of  single 
fiber  and  7x3  fiber  bundle. 


measurably  more  peaks.  There  was  no  statistically  significant  difference  of  the  mean, 
median,  or  mode  between  the  single  fiber  peaks  and  the  fiber  bundle  peaks.  This 
indicates  no  multiple  scatter  events  were  measured  in  the  7x3  fiber  bundle  irradiated 
with  Sr-90.  Similar  results  are  seen  when  the  fibers  are  irradiated  with  Cs-137  and 
Am-241  as  seen  in  Appendix  C. 
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V.  Conclusions 


5.1  Conclusions 

A  MCNP  transport  model  was  developed  with  the  objective  of  analyzing  the 
interaction  of  neutrons  with  fiber  bundles  of  varying  size.  A  2.5  MeV  isotropic  neutron 
source  was  modeled  and  placed  at  different  locations  relative  to  the  fiber  bundles. 
Five  million  source  particles  were  simulated.  A  total  of  42  simulations  were  run,  split 
between  14  simulations  for  each  of  the  four  different  source  locations.  A  gradient  in 
energy  deposition  starts  to  be  seen  starting  at  the  10x10  fiber  bundles  onward.  A 
40x40  bundle  has  6.7±0.2%  intrinsic  efficiency  for  double  scatter  events.  Over  40%  of 
the  double  scatter  events  occur  with  a  time  separtion  of  1  ns  or  greater  in  the  40x40 
bundle. 

In  order  to  measure  only  the  pulses  arising  from  a  PMT’s  dark  current,  an  ab¬ 
solutely  dark  environment  needs  to  be  created  around  the  PMT’s  window.  Three 
methods  at  creating  a  dark  environment  around  the  PMTs  were  studied.  The  lowest 
dark  count  rate  was  achieved  with  the  lens  cap  /  aluminum  can  and  collar  system  with 
a  dark  count  rate  of  15  per  second.  The  other  two  methods  produced  dark  count  rates 
that  were  26x  (for  the  PMT  in  the  large  dark  box)  and  16x  greater  (in  the  aluminum 
dark  box).  Since  a  noisy  environment  is  expected  and  the  scintillation  pulses  may 
be  as  weak  as  dark  count  pulses,  a  method  to  distinguish  between  these  two  types  of 
pulses  is  required.  A  bimodal  distribution  of  decay  constants  was  discovered.  These 
two  distributions  are  statistically  different  from  each  other.  Boundary  values  that 
delineate  dark  pulses  from  scintillation  pulses  can  be  determined,  once  the  distribu¬ 
tion  of  decay  constants  is  known.  A  check  for  these  values  can  be  inserted  directly 
into  the  digital  data  acquisition  code  for  real-time  discrimination  of  dark  events  from 
scintillation  events  using  a  single  PMT. 
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Two  kinds  of  fiber  bundle  were  constructed  for  during  this  effort.  The  first  fiber 
bundle  was  abandoned  for  three  reasons.  It  was  difficult  to  cleanly  apply  the  epoxy, 
the  method  of  cleaving  the  fiber  left  jagged  fiber  faces,  and  there  was  a  slight  bow 
over  the  length  of  the  fiber  bundle.  The  second  type  of  fiber  bundle  created  was  a 
7x3  array  with  no  air  gap  between  the  fibers.  This  fiber  bundle  was  inherently  more 
rigid  because  of  how  the  fibers  are  epoxied  to  each  other,  there  is  no  bowing,  and  the 
fiber  bundle  faces  are  more  reliably  cleaved  and  polished  when  compared  to  the  first 
fiber  bundle  type. 

Measuring  a  scintillation  light  pulse  simultaneously  from  both  ends  of  a  scintillat¬ 
ing  fiber  can  potentially  be  used  to  determine  where  the  scintillation  event  took  place 
along  the  length  of  the  fiber.  This  kind  of  dual-ended  readout  is  also  useful  in  dis¬ 
criminating  scintillation  signals  from  a  noisy  background.  Before  any  data  collection 
using  two  PMTs  took  place,  they  were  gain  matched.  Data  from  dual-ended  readouts 
can  be  analyzed  to  discriminate  out  scintillator  events  from  background  noise.  An 
algorithm  was  developed  to  determine  the  distribution  of  mean  voltages  measured 
from  PMT  A  whenever  no  pulse  from  PMT  B  rose  3  standard  deviations  above  the 
noise  floor.  As  expected,  pulses  measured  from  the  7x3  bundle  are  larger  than  those 
measured  from  the  single  fiber.  This  is  because  the  interaction  medium  is  larger  in 
the  fiber  bundle;  the  light  produced  by  multiple  fibers  should  be  greater  than  that 
produced  by  a  single  fiber.  No  significant  difference  is  measured  between  the  multiple 
peaks  seen  in  scintillation  pulses  of  single  fibers  and  fiber  bundles.  Multiple  scatter 
events  were  not  detected  using  a  7x3  fiber  bundle.  Construction  of  a  40x40  fiber  bun¬ 
dle  is  recommended  for  future  research  due  to  its  intrinsic  efficiency  and  measurable 
double-scatter  event  probability. 
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Appendix  A.  MCNP  Particle  Track  Results 


Figure  40.  MCNP  simulated  energy  deposited  per  fiber  in  a  lxl,  2x2,  3x3,  4x4,  5x5, 
6x6,  7x7,  8x8,  9x9,  10x10,  15x15,  20x20,  40x40,  and  100x100  fiber  bundle  (left  to 
right,  top  to  bottom).  Lighter  colors  represent  more  energy  deposited  than  darker 
colors.  Note  that  the  energy  scale  is  identical  to  that  of  Figure  26.  The  location  of 
the  isotropic  neutron  source  is  defined  by  Part  D  of  Figure  8.  A  noticeable  gradient  in 
energy  deposition  can  be  seen  starting  from  the  15x15  bundle  onward. 
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Figure  41.  MCNP  simulated  energy  deposited  per  fiber  in  a  lxl,  2x2,  3x3,  4x4,  5x5, 
6x6,  7x7,  8x8,  9x9,  10x10,  15x15,  20x20,  40x40,  and  100x100  fiber  bundle  (left  to 
right,  top  to  bottom).  Lighter  colors  represent  more  energy  deposited  than  darker 
colors.  Note  that  the  energy  scale  is  identical  to  that  of  Figure  26.  The  location  of 
the  isotropic  neutron  source  is  defined  by  Part  G  of  Figure  8.  A  noticeable  gradient 
in  energy  deposition  can  be  seen  starting  from  the  20x20  bundle  onward.  Notice  that 
more  energy  is  deposited  deeper  inside  the  bundle  compared  to  Figure  26. 
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Figure  42.  MCNP  simulated  energy  deposited  per  fiber  in  a  lxl,  2x2,  3x3,  4x4,  5x5, 
6x6,  7x7,  8x8,  9x9,  10x10,  15x15,  20x20,  40x40,  and  100x100  fiber  bundle  (left  to 
right,  top  to  bottom).  Lighter  colors  represent  more  energy  deposited  than  darker 
colors.  Note  that  the  energy  scale  is  identical  to  that  of  Figure  26.  The  location  of 
the  isotropic  neutron  source  is  defined  by  Part  H  of  Figure  8.  A  noticeable  gradient 
in  energy  deposition  can  be  seen  starting  from  the  40x40  bundle  onward.  Notice  that 
less  energy  is  deposited  deeper  inside  the  bundle  compared  to  Figure  26. 
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Appendix  B.  Representative  Fiber  Scintillation  Pulses 


Figure  43.  Representative  waveform  measured  from  the  single  fiber  irradiated  with 
Am-241. 


Figure  44.  Representative  waveform  measured  from  the  7x3  fiber  bundle  irradiated 
with  Am-241. 
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Figure  45.  Representative  waveform  measured  from  the  single  fiber  irradiated  with 
Cs-137. 


Figure  46.  Representative  waveform  measured  from  the  7x3  fiber  bundle  irradiated 
with  Cs-137. 
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Appendix  C.  Scintillation  Peak  Distributions 


Figure  47.  Distributions  of  peak  separation  times  measured  during  Cs-137  irradiation 
of  single  fiber  and  7x3  fiber  bundle. 
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Figure  48.  Distributions  of  number  of  peaks  measured  during  Cs-137  irradiation  of 
single  fiber  and  7x3  fiber  bundle. 
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Figure  49.  Distributions  of  number  of  peaks  measured  during  Am-241  irradiation  of 
single  fiber  and  7x3  fiber  bundle. 
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Figure  50.  Distributions  of  peak  separation  times  measured  during  Am-241  irradiation 
of  single  fiber  and  7x3  fiber  bundle. 
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Appendix  D.  Dual  Ended  Readout  Oscilloscope 


1  11  Initialize  Scope 
clc 

clear  all 
close  all 

6  111  Connection  Variables 

visa.brand  =  ’tek’;  1  set  to  ’tek’  for  textronix  scopes 
visa_address  =  ’ TCPIP  :  :  169 . 254 . 235 . 51 :  :  INSTR  ’  ;  1  use  oscilloscope 
ip  address 

buffer  =  4000  *  1024;  7.2000  KiB 

11  111  Open  Instrument 

f  close  (’  all  ’)  ;  1  close  any  open  files  -  useful  if  this  script  was 

interrupted  and  you  want  to  use  it  again 
dpo  =  visa (visa_brand ,  visa_address  ,  ’ InputBuf f er  ’  ,  buffer,  ’... 

OutputBuf f er ’ ,  buffer);  1  create  visa  object,  no  instrument  ... 
driver  required ! 

fopen(dpo);  1  open  the  instrument  so  that  commands  can  be  written 
to  it 

16  f write (dpo ,  ’ *RST 1 ) ;  1  Reset  scope  to  factory  defaults  -  reset  any 
inadvertant  scope  setting  changes, 
f printf (’ Scope  reset  to  factory  defaults  An1)  ; 

II  Configure  Scope 

III  Configuration  Variables 
fprintf ( ’Configuring  scope  .  .  .  .  ’)  ; 

21  frames  =  1000;  1  Number  of  waveforms  to  capture  in  FastFrame  mode 
record  =  1000;  1  Number  of  points  in  each  individual  waveform 
ChlVPos  =  4;  1  from  -5.0  to  5.0  (the  divisions  on  vertical  axis) 
ChlHPos  =  15;  1  from  0  to  100  (divisions  across  the  horizontal  ... 

axis  ) 

Ch2VPos  =  4;  1  from  -5.0  to  5.0  (the  divisions  on  vertical  axis) 

26  Ch2HPos  =  30;  1  from  0  to  100  (divisions  across  the  horizontal  ... 
axis  ) 

ChlScale  =  0.020;  1  [V] 

Ch2Scale  =  0.020;  1  [V] 

ChlOf f  set  =  0.0;  1  [V] 

Ch20f f  set  =  0.0546;  1  [V] 

31  TimeWindow  =  5.0E-9;  ‘/.width  of  time  axis 

TimeStep  =  TimeWindow/record ;  1  Returns  seconds  per  point 

timearray  =  ( 1 : 1 : record )*( TimeWindow / record )* 1 E9 ;  1  Array  of  time 

f write (dpo ,’ ACQ : STATE  O’);  1  Turn  off  acquisition  before  we  ... 

change  scope  settings 

36  f  write  (dpo  ,’ HEAD  O’);  7.  Turn  off  header  information 

f  write  (dpo  ,’ ACQUIRE  :  MODE  SAMPLE’);  7.Sets  sampling  mode  to  Sample  . 
mode 

f write (dpo ,’ ACQUIRE : SAMPLINGMODE  RT ’ ) ;  1  Set  Real  Time  sampling  .. 
mode 
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f  write  (dpo  ,  ’  HOR  :  MODE  MANUAL’);  7»  Enable  manual  adjustment  of  ... 
horizontal  settings 

fprintf  (dpo  HOR :  MODE  RECORD  7,i  ’  ,  record);  7.  Set  number  of  points 
on  horizontal  axis 

41  f  write  (dpo  HOR  :  FAST  :  STATE  1’);  7«  Turn  on  FastFrame  mode 

fprintf  (dpo  HOR  :  FAST  :  COUNT  %i  ’  ,  frames);  7.  Set  number  of  ... 

waveforms  to  capture  per  FastFrame  acquisition 
fprintf  (dpo  DATA  :  START  O’);  7.  Sets  the  beginning  of  the  ... 

individual  waveform  data  to  transfer  (use  0) 
fprintf  (dpo  ,’  DATA  :  STOP  7.i  ’  ,  record);  7«  Sets  the  end  of  the  ... 

individual  waveform  data  to  transfer 
fwrite(dpo,  ’DATA:ENCDG  SRIBINARY’);  7«  Set  data  encoding  to  least 
significant  byte  first,  signed  integer 

46 

fprintf  (dpo  ,’ WFMInpre  :  ENCDG  BINARY’);  7«  Sets  incoming  data  to  ... 
binary  format 

fprintf  (dpo  ,’ WFMInpre  :  BN_FMT  RI  ’  )  ;  7»  Set  signed  integer 
fprintf  (dpo  ,’ WFMInpre  :  BYT_0R  LSB  ’  )  ;  7»Transfer  Least  Significant  .. 
Bit  first 

fprintf  (dpo  ,’  WFMInpre  :  BYT_NR  l’);7«Set  1  byte  format  =  int8 

51 

fprintf  (dpo  ,’ WFMOutpre  :  BYT_NR  1’);  7»Set  1  byte  format  =  intl6 
fprintf  (dpo  ,’ WFMOutpre  :  BN_FMT  RI  ’  )  ;  7»Set  signed  integer 
fprintf  (dpo  ,’ WFMOutpre  :  BYT_0R  LSB’);  7«Transfer  Least  Significant  . 
Bit  first 

56  f  write  (dpo  ,’ Chi  :  Termination  50’)  7«  Sets  input  impedance  in  Ohms  (. 
use  50) 

f  write  (dpo  ,’ Ch2  :  Termination  50’)  7«  Sets  input  impedance  in  Ohms  (. 
use  50) 

f write  (dpo  ,’ Select  :  Chi  On’);  7»  Turns  input  for  channel  1  on 
f write  (dpo  ,’ Select  :  Ch2  On’);  7»  Turns  input  for  channel  2  on 
f  write  (dpo  ,’ DATA  :  SOURCE  Chi,  Ch2  ’  )  ;  7»  Sets  data  source  to  Channel 
1  and  2 

61  fprintf  (dpo  ,’  Horizontal  :  Position  7»i’,  ChlHPos);  7«  Places  signal  at 
designated  horizontal  location 

fprintf  (dpo  ,’ chi  :  Scale  7,i  ’  ,  ChlScale);  l  Sets  vertical  scale  for 
channel  1 

fprintf  (dpo  ,’ ch2  :  Scale  %i  ’  ,  Ch2Scale);  7,  Sets  vertical  scale  for 
channel  2 

fprintf  (dpo  ,’ CHI  :  Position  7«i’,  ChlYPos);  7o  Aligns  channel  1 

fprintf  (dpo  ,’ CH2  :  Position  7«i’,  Ch2VPos);  7»  Aligns  channel  2 

66  fprintf  (dpo  ,’  Chi  :  Of  f  set  7oi’,  ChlOffset);  7.  Sets  channel  1  offset 

fprintf  (dpo  ,’ Ch2  :  Of  f  set  7oi’,  Ch20ffset);  7.  Sets  channel  2  offset 

fprintf  (dpo  ,’ Horizontal  :  Scale  7oi  ’  ,  TimeWindow);  7»  Sets  time  scale 
fprintf  (dpo  ,’ Horizontal  :  AcqLength  7»i  ’  ,  record);  7o  Sets  time  ... 
resolut ion 

71  query ( dpo,  ’WFMOutpre: YMult  ? ’ )  ; 

ChlYMult  =  st  r2double  (  query  (  dpo  ,’ WFMOutpre  :  YMult  ?’))  ;  7.  Returns  . 

vertical  scale  factor 

ChlYoff  =  -ChlVPos  *ChlScale  ;  7«  Calculate  offset  from  zero 


Ch2YMult  =  s tr2double ( query ( dpo ,  ’ WFMOutpre : YMult ? 1 )) ;  1  Returns 

vertical  scale  factor 

Ch2Yoff  =  -Ch2VPos *Ch2Scale ;  7«  Calculate  offset  from  zero 

76 

7. 7.7. 1  "/.Edge  Trigger  Sett  ing  7,7, 7.7.7. 

f write (dpo ,’ Trigger : A : Type  Edge’);  7.  Edge  mode  trigger 
f  write  (dpo  ,  ’  Trigger  :  A  :  Edge  :  Slope  Fall’);  7.  Falling  slope 
81  f  write  (dpo  ,  ’  Trigger  :  A  :  Level  -0.010’);  7.  Set  trigger  voltage  level 
fprintf  (dpo  Trigger  :  A  :  Mode  Normal’);  7.  Set  trigger  to  Normal,  .. 
apparently  "Auto"  is  bad 


HmPulse  Width  Trigger  Setting 7» 7« 7. 1 
86  /{ 

f  write  (dpo  ,’ Trigger  :  A  :  Mode  Normal’);  7.  Set  trigger  to  Normal,  ... 
apparently  "Auto"  is  bad 

f  write  (dpo  ,’ Trigger  :  A  :  Type  Pulse’);  7»  Pulse  mode  trigger 
f  write  (dpo  ,’ Trigger  :  A  :  Pulse  :  Class  Width’);  7.  Set  to  Pulse  Width  .. 
mode 

f  write  (dpo  ,’ Trigger  :  A  :  Pulse  :  Width  :  HighLimit  10.3E-9’);  7.  Set  upper 
pulse  limit  [s] 

91  f  write  (dpo  ,’ Trigger  :  A  :  Pulse  :  Width  :  LowLimit  6.2E-9’);  7.  Set  lower  . 
pulse  limit  [s] 

f  write  (dpo  ,’ Trigger  :  A  :  Pulse  :  Width  :  Polarity  Negative’);  7.  Trigger  . 
off  of  negative  polarity  pulses 

f  write  (dpo  ,’ Trigger  :  A  :  Pulse  :  Width  :  When  Within’);  7»  Trigger  when  .. 

the  time  window  occurs  within  the  pulse 
f  write  (dpo  ,’ Trigger  :  A  :  Level  -0.04’);  7.  Set  trigger  voltage  level; 

7.} 
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fprintf ( ’Done  !\n’)  ; 
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7.7.  ACQUIRE  DATA 

sets  =  1000;  1  How  many  datasets  do  you  want  to  record? 

tic  7.  Start  run  time  clock 
Chi  =  zeros  (1000 , 1000)  ; 

Ch2  =  zeros  (1000 , 1000)  ; 
for  i=l : sets 

fprintf (dpo ,’ Acquire : StopAf ter  Sequence’);  1  Stop  acquiring  .. 

after  FastFrame  acquisition 
f write (dpo ,’ ACQ : STATE  1’);  1  Turn  on  acquisition 

fprintf(’Acquiring  data  set  #7.  i . ’,  i); 

dpo. Timeout  =  90;  1  Adjust  MATLAB  timeout  so  it  is  longer  than 
the  acquisition  duration 

query  ( dpo  ,  ’*  opc  ?’)  ;  7.  Will  return  "1"  when  the  acquistion  is  . 

complete  -  used  to  sync  curve  command 
dpo. Timeout  =  10;  1  Decrease  timeout  time  now  that  the  ... 

acquisition  is  done 
fprintf (’Done !\n’) ; 
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f printf ( ’  Asking  for  waveforms....’); 

fwrite(dpo,  ’curve?1);  °i  Ask  for  waveforms  -  waveform  preamble 
and  binary  data  will  be  stored  in  memory 

116  f printf ( ’ Done ! \n ’ ) ; 

fprintf ( ’  Reading  waveforms  ....’); 

f  read  (dpo  ,  1)  ;  7«  Read  and  throw  away  first  character  in  ... 

preamble  -  should  be  a  "#" 

bytes  =  st r2double  (  char  ((  f  read  ( dpo  ,  1 ))))  ;  7«  Read  and  store  the 
number  of  bytes  in  memory 

121  fread(dpo,  bytes);  70  Read  and  throw  away  the  appropriate  ... 

number  of  characters  remaining  in  the  preamble 

header  =  (1+1+bytes); 

rawl  =  fread(dpo,  record*f  rames  ,  ’  int8  ’  )  ;  7«  Read  and  store  all 

of  the  waveform  data  -  no  separation  between  individual  ... 
wavef  orms 

fread(dpo,  header); 

raw2  =  fread(dpo,  record*f  rames  ,’  int8  ’)  ;  7,  Read  and  store  all 
of  the  waveform  data  -  no  separation  between  individual  ... 
wavef  orms 

126  f printf (’ Done ! \n ’)  ; 

f  read  (dpo  ,  1)  ;  7«  Read  and  throw  away  the  termination  character 

that  Curve?  adds 


start  =  1 ; 
stop  =  1000; 

131  for  j=l:frames 

Chi ( 1 : r e cord  ,  j)  =  rawl ( st art : stop )  ; 
Ch2 ( 1 : r e cord  ,  j)  =  raw2 (start : stop) ; 
start  =  stop+1 ; 
stop  =  start +  record -  1  ; 

136  end 


Chi  =  -  (  Ch  1  YMult  *  Chi  +  Chl  Yof  f  )  ;  7»  Apply  scale  and  offset  to  ... 

match  the  physical  voltage  values 
Ch2  =  -  ( Ch2YMult  *  Ch2  +  Ch2Yof  f )  ;  7»  Apply  scale  and  offset  to  ... 

match  the  physical  voltage  values 
savename  =  [’Chl_’  num2str(i)  ’.mat’];  '/,  name  the  current  ... 
batch  of  waveforms 

141  save  ( savename  ,’ Chi  ’  )  7«  save  the  current  batch  of  waveforms 

savename  =  [  ’  Ch2_  ’  num2str(i)  ’.mat’];  7.  name  the  current  ... 
batch  of  waveforms 

save ( savename ,’ Ch2 ’ )  i  save  the  current  batch  of  waveforms 

end 

f  pr  int  f  (’*********  *  7oi  waveforms  recorded!  =t==t=>t==i==i==t:=(==(==^=(=\n’  ,  sets*, 
frames  *2)  ; 

146  timeElapsed  =  toe  °/  End  run  time  clock 
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Appendix  E.  Particle  Track  Code 


2  MCNP5  Fiber  Bundle  Particle  Track  Parser 

Given  a  set  of  particle  track  files  with  the  filename  structure  .. 
given 

by:  #x#_@.ptr,  where  #  is  the  x  and  y  dimensions  (respectively  of 
the  fiber 

bundle  and  @  is  a  character  corresponding  to  the  data  set  being  .. 
processed 

n 

7  close  all 
clear  all 
clc 

l  this  array  corresponds  to  the  name  of  the  mcnp  ptrac  files  being 
12  */.  processed 

dim  =  [1,2,3,4,5,6,7,8,9,10,15,20,40,100]; 
data.sets  =  [ ’ a ’  ,  ;  b  ’  ,  ’  c  ’  ,  ’  d  ’  ,  ’  e  ’  ,  ’  f  ’  ,  ’ g ’  ,  ’  h  ’  ]  ; 

for  r  =  1 : length (data_sets ) 

17  set  =  data.sets  (r )  ;  °/„  data  set  being  processed  (the  in  the 

introduction) 

for  q  =  1 : length (dim)  %  loop  through  all  files 

fprintf  (’  Processing  :  yoix°/0i_0/oC  .  ptr\n  ’  ,  dim(q),  dim(q),  set); 
nscatter  =  zeros  ( 1  ,  100)  ; 
int  _  1  i  st  =  []; 

22  energy  =  []  ; 

tic;  °/0  start  run  time  clock 
V.  open  file 

fid  =  fopen([set  ’/’  num2str (dim (q) )  ’x’  num2str (dim (q) )  . 

set  ’  . ptr  ’  ] )  ; 
for  i  =  1  :  1 0 

27  fgetl(fid);  ‘/.skip  first  10  lines 

end 

line  =  1;  °/n  current  line  number  being  looked  at,  starts  after 

header 

k  =  1;  %  index  used  to  keep  track  of  energy  deposited 

m  =  1;  '/.  index  used  to  keep  track  of  interaction  events 

32  hits  =  0;  °/n  number  of  histories  that  occur  within  bundle 

°/.{ 

The  while  loop  reads  from  the  file  one  line  at  a  time  until  it 
reaches  the  end  of  the  file.  It  continuously  builds  up  the 
variable  called  "data"  which  is  just  a  list  of  the  individual 
lines 

37  read  in.  When  it  encounters  a  "9000"  in  the  first  column  it  . 

reads 

the  next  line  and  stops  to  perform  some  calculations.  9000 
corresponds  to  a  terminal  event  in  the  particle  track.  The  .. 
calcuations 

at  this  point  create  the  interaction  list  and  track  the  hits 

n 
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while  ( ~  f  eof  ( f  id )  )  °/0  read  from  file  until  you  reach  the  end  of 
the  file 

s  =  fgetl(fid);  7,  read  a  line  from  the  file 

char.line  =  sscanf  (s(:)  ,  ’7,c’)  ;  7.  read  all  characters  . 

in  line 

data  (  line  ,  1  :  length  (  char  _1  ine  )  )  =  char_line;  7,  buildup  list 

of  lines 

if  ( strcmp (data(line,7:ll),  ’  9000’)  ==  1) 

47  s  =  fgetl(fid)  ;  7,  read  another  line  from  the  file 

char_line  =  sscanf  (s  (  :  )  ,  ’7c’);  7,  read  all  characters 

in  line 

data(line  +  l,l:length(char_line))  =  char_line;  7«r  ead  .. 
last  line 

hits  =  hits+1;  7.  increment  hits 

data_size  =  size(data);  7»  calculate  dimensions  of  data 
variable 

52  data_rows  =  data.size  (1)  ;  7«  pull  out  #  rows  from  size 

variable 

7.{ 

the  for  loop  cycles  through  the  current  data  block  ... 
until  it 

finds  a  4000,  corresponding  to  a  collision  event,  it  . 
then 

record  the  cell  number,  energy  deposited  and  the  atom 
the 

57  collsion  took  place  with  into  the  interaction  list 

7.} 

for  i=l : data.rows 

if  st r cmp ( dat a ( i  ,  7 : 1 1 )  ,  ’  4000')  ==  1 

int_list(m,  1)  =  str2double ( data ( i  ,  46 : 5 1 ) )  ; 

62  int_list(m,  2)  =  st r 2double ( dat a ( i +2  ,  28 : 32 ) )  ; 

int_list(m,  3)  =  st r 2double ( dat a ( i  + 1  ,  8 1 : 93) )  . 

-  str2double ( data ( i +3 , 81 : 93) )  ; 
int_list(m,  4)  =  m; 
m  =  m+ 1 ; 

67  end 

end 

p_energy  =  str2double (data (data.rows  ,  81 : 93) )  ; 
collision  =  str2double ( data ( data.rows - 1  ,  69 : 71 ) )  ; 
if  isnan ( collision)  ==  0 
72  if  (collision  ~=  0) 

nscatter ( collision)  =  nscatter ( collision) +1 ; 

end 

end 

if  (p_energy  ~=  2.5) 

77  energy(k)  =  2.5  -  p_energy; 

k  =  k+1; 

end 

clear  data; 
line  =  1 ; 

82  end 
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line 


line  +  1  ; 


end 

save ( [set  ’/’  num2str (dim (q) )  ’ x  ’  num2str (dim (q) )  ’_’  set  ’. 

mat ’  ]  ,  ... 

’int_list’,  ’energy’,  ’nscatter’,  ’hits’) 

87  toe 

end 
end 
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Appendix  F.  MCNP  Input  File  Generation  Code 


1  close  all 
clear  all 
clc 

dim  =  [1,2,3,4,5,6,7,8,9,10,15,20,40,100]; 
set  =  ’g  ’  ; 

6  for  q  =1  :  length,  (dim) 
tic  ; 

number_f  ibers.x  =  dim(q);  7,  number  of  fibers  in  the  x-direction 
number.f  ibers_y  =  dim(q);  70  number  of  fibers  in  the  y-direction 
f  iber.diameter  =  0.0500;  7o  fiber  diameter  [cm] 

11  fiber_length  =  19.0*0.5;  7,  fiber  length  in  the  z-direction  [cm] 

total_fibers  =  number  _f  ibers.x  *  number  _f  ibers  _y  ;  7.  total  number  of 

fibers  in  the  bundle 

spacing  =  f  iber  _di  amet  er  *  1 . 0 1  ;  7«  cells  that  hold  fibers  are  7ol  ... 

bigger  than  the  fibers  themselves 

world_diameter  =  200;  7»  diameter  of  world  sphere  [cm] 

16  world_id  =  t  otal_f  ibers  +  number  _f  ibers.x  +  number_f  ibers_y  +2  +  3  ;  7»  ... 

surface  id  for  world  sphere 

fiber_start_id  =  total_f iber s  +  number _f iber s_x  +  number _f ibers_y  +  2+l ; 
‘/,  id  of  start  plane  for  fibers 

fiber_end_id  =  f  iber_start_id  +  l  ;  7«  id  of  end  plane  for  fibers 

source.pos  =  [-5,  number _f iber s _y * f iber _diamet er +0 . 5 ,  fiber_length 

*0.5];  7,  xyz  position  of  source,  fiber  bundle  is  in  first  ... 

quadrant,  positive  z 

21  number.of .particles  =  5000000;  %  number  of  particles  to  run 


x_planes  =  zeros  (  number  _f  iber  s_x  +  1  ,  1 )  ;  7«  planes  that  form  borders 

between  fiber  cells  in  x-direction 
y_planes  =  zeros  ( number  _f  iber  s_y +  1  ,  1 )  ;  7«  planes  that  form  borders 

between  fiber  cells  in  y-direction 


26  x.planes  (1)  =  -spacing/2;  7«  first  border  plane  in  x-direction 

for  k=2 : number.f ibers_x+l 

x_planes(k)  =  x_planes (k-1) +spacing ; 

end 

y_planes(l)  =  -spacing/2;  "/,  first  border  plane  in  x-direction 
31  for  k=2 : number.f ibers_y +1 

y_planes(k)  =  y_planes (k-1) +spacing ; 

end 


fid  =  f open ( [num2str (number.f ibers.x)  ’x’  num2str ( number _f iber s _y ) 
’  _  ’  set  ’.in’],  ’  w  ’  )  ;  7.  open  file  for  writing 

36  fprintf(fid,  ’fiber  array  test\n’);  7»  title 
for  i  =  1  :  t  ot  al_f  iber  s  7»  write  fiber  cells 

fprintf(fid,  ’7oi  1  -1.05  -7,i  7oi  ~7oi  imp:n  =  l\n’,  i,  i,  ... 
f iber_start_id ,  f iber_end_ id ) ; 

end 
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41  n  =  total_fibers  +  l;  °L  current  cell  index  after  writing  fiber  cells 
for  i  =  1  +  number _f ibers_y  +  total_ fibers : number _f ibers_y  +  .  .  . 
number _f ibers_x+total_f ibers 
for  j=l+total_f ibers: number _f ibers_x+total_f ibers 

fprintf (fid,  ’‘/,i  0  °/ni  '/.i  -°/0i  %i  -%i  '/.i  -°/0i  imp:n=l\n’,  n, 
n-total_f ibers  ,  j,  j+1,  i  +  1,  i+2,  f iber_start_id  ,  ... 

f iber _end_ id ) ;  %  fiber  cells 

n  =  n  +  1;  ‘L  increment  cell  index 

46  end 

end 


fprintf(fid,  ,0/ni  0  #(°/„i  -‘Li.  ‘Li.  ~‘Li-  ‘Li.  ~°L  i)  -'/.i  imp  :  n  =  l\n 1  ,n ,  1+... 
total_f ibers  ,  t ot al _f iber s +  number_f iber s_x +  1 ,  t ot al_f iber s  +  .  .  . 
number _f ibers_x  +  2 ,  2+ number _f iber s_y  +  number _f ibers_x  +  . .  . 

total_f ibers ,  f iber_start_id ,  f iber_end_id ,  world_id);  %  cell 
that  holds  bundle 

51  fprintf (fid,  ’7,i  0  %i  imp  :  n  =  0\n’,  n  +  1,  world_id);  ‘L  world  cell 


fprintf(fid,  ’Xn1); 

fprintf (fid,  ’c  ===  Surface  Cards  ... 

=======================================================\n’ ) ; 

n  =  1 ; 

56  for  y =0 : spac ing :( number _f iber s_y - 1 )* spac ing 

for  x=0: spacing: ( number _f ibers_x -1) ^spacing 

fprintf (fid,  ’°/0i  c/z  %f  °/0f  %f\n’,  n,  x,  y  ,  ( f  iber.diameter . 
/  2  )  )  ; 
n  =  n+ 1 ; 

end 

61  end 


for  i =1 : length ( x_plane s ) 

fprintf(fid,  ,0/0i  px  70f\n’,  n,  x.planes  ( i )  )  ; 
n  =  n+ 1 ; 

66  end 

for  i = 1 : length ( y_plane s ) 

fprintf(fid,  ,0/0i  py  70f\n’,  n,  y_planes  ( i )  )  ; 
n  =  n+ 1 ; 
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end 

fprintf (f id  , 

’  ‘Li.  pz 

0\n  ’  , 

f iber_s 

fprintf  (f id  , 

’  L  i  pz 

7,  f  \  n  ’ 

,  fiber_ 

fprintf (f id  , 

’  ‘L  i  so 

7.f\n> 

,  world. 

sphere 

fprintf (f id  , 

1 \n  ’  )  ; 

fprintf (f id  , 

’  c  === 

Materials  . .  . 

=  \n  ’  ) 


fprintf(fid,  ;ml  1001  0.50  6012  O.SOXn1); 

fprintf(fid,  ’c  ======  Source  ... 

===========================================================\n0 
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fprintf(fid,  ’  sdef  erg  =  2.5  par  =  l  pos  =  ‘Li  %f  7«f\n’,  source.pos  (1)  ,  .  .  . 

source.pos (2)  ,  source.pos  (3) )  ; 
fprintf(fid,  ’c  ======  Histories  ... 

=================================================== =====\n. ’ ) ; 

81  fprintf(fid,  ’nps  7of\n’,  number.of .particles)  ; 

fprintf (fid,  ’c  ======  Cutoffs,  physics,  and  variance  reduction... 

===========================\n’); 

fprintf (fid,  ’mode  n\n’); 

fprintf (fid,  ’ptrac  write  =  all  file  =  asc  max=7„f  &\n’,  ... 

number.of .particles *20) ; 

fprintf(fid,  ’  filter  =  l,"/0i,icl\n’,total_fibers); 

86  '/,  fprintf (fid,  ’  cell=’); 

7,  for  i  =  l:total_fibers-l 
7.  if  (mod(i,10)  ~  =  0) 

7.  fprintf  (f  id  ,  ’7.i,’,  i); 

7.  else 

91  7«  fprintf  (fid,  ’  &\n  7oi  ,  ’  ,  i); 

7c  end 

7c  end 

7c  fprintf  (fid,  ’  7«  i  \  n  ’  ,  total.f  ibers)  ; 

96  f close  all ; 

toe 

end 
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Appendix  G.  Find  Maxima  Code 


1  function  [xmax  ,  imax  ,  xmin  ,  imin]  =  extrema(x) 

"/.EXTREMA  Gets  the  global  extrema  points  from  a  time  series. 

/  [XMAX , IMAX , XMIN , IMIN]  =  EXTREMA (X)  returns  the  global  minima  . 
and  maxima 

7.  points  of  the  vector  X  ignoring  NaN’s,  where 
7.  XMAX  -  maxima  points  in  descending  order 
6  7.  IMAX  -  indexes  of  the  XMAX 

7.  XMIN  -  minima  points  in  descending  order 
i  IMIN  -  indexes  of  the  XMIN 

7. 

7.  DEFINITION  (from  http://en.wikipedia.org/wiki/... 
Maxima_and_minima ) : 

11  7.  In  mathematics,  maxima  and  minima,  also  known  as  extrema,  are 
points  in 

7,  the  domain  of  a  function  at  which  the  function  takes  a  largest 
value 

7,  (maximum)  or  smallest  value  (minimum),  either  within  a  given 
7.  neighbourhood  (local  extrema)  or  on  the  function  domain  in  its 
ent irety 

7.  (global  extrema). 

16  7. 

7.  Example  : 

7.  x  =  2*pi*linspace(-l,l); 

7.  y  =  cos(x)  -  0.5  +  0 . 5*  rand  (  size  (x)  )  ;  y(40:45)  =  1.85;  y... 

(50 : 53) =NaN ; 

7.  [ymax  ,  imax  ,  ymin  ,  imin]  =  extrema  (y)  ; 

21  7.  plot (x , y , x ( imax )  , ymax ,  ’  g  .  5  , x ( imin )  , ymin  ,  ’  r  .  5  ) 

1 

7.  See  also  EXTREMA2  ,  MAX,  MIN 
7.  Written  by 

26  7.  Lie.  on  Physics  Carlos  Adrin  Vargas  Aguilera 
7.  Physical  Oceanography  MS  candidate 
7.  UNIVERSIDAD  DE  GUADALAJARA 
/  Mexico ,  2004 

1 

31  7.  nubeobscura@hotmail.com 


http : // www . mathworks . com/m at  lab central/ f ileexchange 
12275 

2006-09-14 

English  translation  from  Spanish. 

Accept  NaN  ’  s  . 

Change  name  to  MAXIMA,  and  definition  added. 

41  xmax  =  []  ; 

imax  =  []  ; 

xmin  =  []  ; 

imin  =  []  ; 


7.  File  ID 
7.  Submited  at 
36  7.  2006-11-11 
7.  2006-11-17 
2007-04-09 
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46 

51 

56 

61 

66 

71 

76 

81 

86 

91 


%  Vector  input? 

Nt  =  numel (x)  ; 
if  Nt  ~=  length(x) 
error ( ' Entry  must  be  a  vector.’) 
end 


7„  NaN  ’  s  : 

inan  =  f ind ( isnan  (x) )  ; 
indx  =  1 : Nt ; 
if  ~ isempty ( inan) 
indx  (  inan  )  =  []  ; 

x  (  inan)  =  []  ; 

Nt  =  length (x) ; 
end 


“/«  Difference  between  subsequent  elements: 
dx  =  dif f  (x)  ; 


1  Is  an  horizontal  line? 
if  ~  any ( dx ) 
return 
end 


1  Flat  peaks?  Put  the  middle  element: 


a  =  f ind ( dx  ~  =  0 )  ; 
lm  =  f  ind  (  dif  f  (  a)  ~  =  1 )  +  1; 

d  =  a(lm)  -  a(lm-l) ; 
peak 

a(lm)  =  a(lm)  -  floor(d/2) 
a ( end+1 )  =  Nt ; 


“/,  Indexes  where  x  changes 

“/„  Indexes  where  a  do  not  changes 

%  Number  of  elements  in  the  flat 

%  Save  middle  elements 


“/.  Peaks? 
xa 
b 


xb 

imax 

imin 

imax 

imin 


=  X  (  a)  ; 

1  S  e  r  i  e 

without  flat  peaks 

(dif f (xa)  >  0)  ; 

1  1 

=  > 

positive  slopes  (minima 

begin ) 

1  0 

=  > 

negative  slopes  (maxima 

begin ) 

=  dif f (b)  ; 

1  -1 

=  > 

maxima  indexes  (but  one) 

1  +1 

=  > 

minima  indexes  (but  one) 

=  f ind ( xb  ==  - 1 ) 

+  1 ;  °/« 

maxima  indexes 

=  f ind ( xb  ==  + 1 ) 

+  1 ;  °/» 

minima  indexes 

a ( imax ) 
a ( imin ) 


nmaxi  =  length ( imax ) ; 
nmini  =  length ( imin)  ; 


*/.  Maximum  or  minumim  on  a  flat  peak  at  the  ends? 
if  (nmaxi==0)  &&  (nmini==0) 
if  x(l)  >  x(Nt) 
xmax  =  x  (  1 )  ; 
imax  =  indx ( 1 ) ; 
xmin  =  x ( Nt )  ; 
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96 

101 

106 

111 

116 

121 

126 

131 

136 

141 

146 


imin  =  indx ( Nt ) ; 
elseif  x(l)  <  x(Nt) 
xmax  =  x ( Nt )  ; 
imax  =  indx ( Nt ) ; 
xmin  =  x  (  1 )  ; 
imin  =  indx  (  1 )  ; 
end 

return 

end 

*/.  Maximum  or  minumim  at  the  ends? 
if  (nmaxi==0) 
imax (1:2)  =  [1  Nt] ; 

elseif  (nmini==0) 
imin (1:2)  =  [1  Nt] ; 
else 

if  imax(l)  <  imin(l) 
imin (2 : nmini + 1 )  =  imin; 

imin  (  1 )  =  1 ; 

else 

imax (2 : nmaxi + 1 )  =  imax; 

imax ( 1 )  =  1 ; 

end 

if  imax(end)  >  imin(end) 
imin(end+l)  =  Nt ; 
else 

imax(end+l)  =  Nt ; 
end 
end 

xmax  =  x ( imax ) ; 
xmin  =  x ( imin )  ; 

°i  NaN  ’  s  : 

if  ~ isempty ( inan) 
imax  =  indx(imax); 
imin  =  indx(imin); 
end 

“/o  Same  size  as  x: 

imax  =  r eshape ( imax , size ( xmax )) ; 
imin  =  reshape ( imin , size (xmin) ) ; 

“/.  Descending  order: 

[temp , inmax]  =  sort (-xmax);  clear  temp 
xmax  =  xmax (inmax); 
imax  =  imax(inmax); 

[xmin , inmin]  =  sort (xmin); 
imin  =  imin(inmin); 

7o  Carlos  Adrin  Vargas  Aguilera,  nubeobscura@hotmail.com 
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